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Abstract 

This paper presents DELPHI measurements and interpretations of cross- 
sections, forward-backward asymmetries, and angular distributions, for the 
e’''e“ —ff process for centre-of-mass energies above the Z resonance, from 
a/s ~ 130 — 207 GeV at the LEP collider. 

The measurements are consistent with the predictions of the Standard Model 
and are used to study a variety of models including the S-Matrix ansatz for 
e’''e“ —> ff scattering and several models which include physics beyond the Stan¬ 
dard Model: the exchange of Z bosons, contact interactions between fermions, 
the exchange of gravitons in large extra dimensions and the exchange of v in 
R-parity violating supersymmetry. 
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1 Introduction 

This paper presents measurements and interpretations of cross-sections, cr, forward- 
backward asymmetries, Aps, and angular distributions, da/dcos^, for e+e“ —ff pro¬ 
cesses for centre-of-mass energies above the Z resonance, as measured in the DELPHI 
experiment at the LEP collider. Measurements of flavour tagged qq samples will be 
included in an additional publication. 

For the hrst part of the LEP programme, LEP I, e’''e“ collisions were provided at 
centre-of-mass energies close to the Z resonance. Measurements of the process e+e“ —ff 
were used to determine properties of the Z and electroweak parameters of the Standard 
Model [1,2]. In 1995 the operation of LEP moved into the LEP II programme. The 
collision energy was raised to signihcantly above the Z resonance, and a total of approx¬ 
imately 0.7 fb“^ of integrated luminosity was delivered to the DELPHI experiment at 
energies ranging from 130 GeV to a maximum collision energy of 209 GeV during LEP 
H operations. A breakdown of the centre-of-mass energies, and integrated luminosities is 
given in Table 1. Measurements of the process e’''e“ —ff from LEP H are less sensitive 
to the electroweak parameters of the Standard Model. Nevertheless, taken together, the 
measurements constitute a test of the Standard Model at the 0{1%) level, at the highest 
e’''e“ collision energies to date. Furthermore, the e+e“ —> ff measurements at LEP H 
are predicted to be more sensitive to a variety of models of physics beyond the Standard 
Model than the LEP I measurements. Having determined many of the parameters of the 
Standard Model, largely from LEP I data, studies of the process e+e“ —ff at LEP H 
are, therefore, a suitable place to look for physics beyond the Standard Model. 

Results from the analysis of data at centre-of-mass energies from 130-189 GeV have 
already been published [3,4]. Results at energies from 192-207 GeV are published here for 
the hrst time. Some previously published results at energies of 183 and 189 GeV [4] have 
been reanalysed and are presented again here. For completeness we present again the 
results from 130-172 GeV data which have not been reanalysed. Results on ff production 
from the other LEP experiments can be found in [5-7]. 

The operation of LEP during the LEP H programme is discussed in Section 2. The 
DELPHI detector itself is described in Section 3, and the measurement of the integrated 
luminosity is described in Section 4. Features of ff production at LEP H, details of the 
event simulation and theoretical predictions are given in Section 5. Sections 5.1 to 5.4 
cover the analysis of the individual channels e+e“ —> e’''e“, e^e“ —, e’''e“ — 
and e+e“ —qq. In each section the criteria for selecting events are described and the 
methods for evaluating the efficiency and backgrounds are discussed. Results of the 
individual measurements are provided and the principal sources of systematic error are 
described. The sets of measurements are summarised in Section 5.5, where the results are 
compared to the predictions of the Standard Model. Further interpretations of the results 
are presented in Section 6. In particular the results, together with LEP I data, are inter¬ 
preted with the S-Matrix formalism in Section 6.1, and also in a variety of models which 
include explicit forms of physics beyond the Standard Model; models with Z bosons 
in Section 6.2, contact interactions (Section 6.3), models which include the exchange of 
gravitons in large extra dimensions (Section 6.4) and models which consider possible s 
or t channel sneutrino Ui exchange in R-parity violating supersymmetry (Section 6.5). In 
each case, a resume of the model, the method of comparing predictions of the model to 
the data and the results of the interpretation are provided. Gonclusions drawn from the 
DELPHI analyses of the e^e“ —ff processes at LEP H are given in Section 7. 



Year 

Nominal 

Energy 

(GeV) 

Mean 

Energy 

(GeV) 

Integrated 

Luminosity 

(pb"^) 

1995 

130 

130.3 

3 


136 

136.3 

3 

1996 

161 

161.5 

10 


172 

172.4 

10 

1997 

183 

182.7 

53 

1998 

189 

188.6 

155 

1999 

192 

191.8 

25 


196 

195.7 

76 


200 

199.7 

83 


202 

201.8 

40 

2000 

205 

204.9 

82 


207 

206.6 

135 

Total 


195.6 

675 


Table 1: Nominal centre-of-mass energies with approximate centre-of-mass energy and 
integrated luminosity collected by DELPHI. The values used in the final analyses differ 
due to data collection efficiencies and selection criteria. 
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2 LEP 

The LEP collider was upgraded from its original configuration, used for running at 
energies around the Z resonance, by the addition of superconducting RF cavities. This 
then allowed operations at energies well above the Z resonance. In the years 1995 to 1999 
LEP delivered e^e“ collisions at one or more discrete energies, each LEP £11 corresponding 
to a particular energy. In 2000 the luminosity at any given time was delivered at the 
maximnm energy available from the LEP RF system, within a margin of safety which 
allowed for the trip of either two, one or no RF nnits before the beam was lost. In a 
given fill the energy would usually be increased from the limit set by two RF trips, to 
the limit set by one trip and eventually no RF trips. The Inminosity was, therefore, 
delivered more or less continnously over a range of energies. For analysis pnrposes the 
data were gronped into two energy points: data taken at centre-of-mass energies between 
202.5 and 205.5 GeV and data taken at energies above 205.5 GeV. Data collected dnring 
the time in which the energy was increasing are not analysed here, which represents a 
loss of approximately 1% of the delivered luminosity. 

As well as providing collisions at energies above the Z resonance, LEP also ran at a 
centre-of-mass energy close to the pole of the Z resonance in each year. The data gathered 
at this energy were used by DELPHI for detector calibration pnrposes. Typically 2.5 pb~^ 
were delivered at the start of the year in 1996-2000, with additional luminosity delivered 
on several other occasions when reqnested by the experiments. The data collected in 
1995 at centre-of-mass energies of 130 and 136 GeV followed a long rnn at energies close 
to the Z resonance. 

The energy of the e’''e“ collisions was determined by the LEP energy gronp. Dnring 
LEP I this was based primarily on the resonant depolarisation techniqne [8], where the 
energy was determined with very high precision at the end of a large number of fills - 
making the measurements at close to the actual collision energies. The measurements 
were used to normalise a model of the RF system, from which the energies during actual 
collisions were determined. At LEP II the energy was again calibrated using the resonant 
depolarisation techniqne, bnt it was not possible to obtain polarisation at the actnal colli¬ 
sion energies, so the RF model had to be used to extrapolate from the calibration energies 
(np to 60 GeV per beam) to the collision energies. The accnracy of this extrapolation 
was checked using a nnmber of techniques: 

• dedicated LEP runs, at collision energies, in which the energy of the electron and 
positron beams were determined from the deflection of the beams in a magnetic 
spectrometer whose magnetic field was known to high precision; 

• measnrements of the synchrotron tune as a fnnction of RF voltage; 

• measnrements of the total magnetic field seen by the beams using the finx loop 
system of LEP. 

The precision that was obtained was ±(20 — 40) MeV on the centre-of-mass energy [9]. 
The instantaneous difference in energy between the electron and positron beams was less 
than ±100 MeV. 

The beam energy spread at LEP II is larger than at LEP I, dne to the increased 
synchrotron energy loss at the higher beam energies. At LEP I the beam energy spread 
was ~ 40 MeV, at LEP II this increases to ~ 180 MeV at the highest centre-of-mass 
energies. At LEP I the beam energy spread had to be taken into acconnt in the estimation 
of the total cross-section. For example, at the peak of the Z resonance, the expected 
cross-section was significantly lower than wonld have been observed for a monochromatic 
beam. At LEP II the corresponding correction is not significant for e^e“ —ff since 
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the cross-sections are nearly linear over a small range of energies. However, for certain 
measurements at LEP II it is necessary to account for the beam energy spread for the 
determination of the event kinematics. 


3 The DELPHI Detector 

A detailed description of the DELPHI apparatus as used at LEP I and its performance 
can be found in refs. [10,11]. For the present analysis the following parts of the detector 
are relevant: 

• for the measurement of charged particles the Microvertex Detector (VD), the Inner 
Detector (ID), the Time Projection Chamber (TPC), the Outer Detector(OD) and 
the Forward Chambers A and B (FCA and FCB). For the running from 1995 onwards 
a lengthened Inner Detector was installed. The polar angle^ coverage was thereby 
extended from 23° < 6 < 157° to 15° < 9 < 165° with a corresponding increase 
in forward tracking efficiency. For a period in 2000 when part of the TPC was not 
operational (see later), space points from the Barrel RICH detector were included 
in the track £t; 

• for the measurement of electromagnetic energy the High-density Projection Chamber 
(HPC) and the Forward Electromagnetic Calorimeter (FEMC); these detectors were 
also used for identifying minimum ionising particles; 

• for the measurement of the hadronic energy and muon identification the Hadron 
Calorimeter (HCAL), which covered both the barrel and endcap regions; 

• for muon identification the barrel (MUB), the endcap (MUF), and the surround 
muon chambers (SMC), which completed the polar angle coverage between barrel 
and endcap; 

• for the trigger, besides the detectors mentioned above, the barrel Time of Flight 
counters (TOF), the endcap scintillators (HOF) and a scintillator layer embedded 
in the HPC; 

• for the measurement of luminosity the Small Angle Tile Calorimeter (STIC). 

The DELPHI detector was upgraded for LEP H data taking. Changes were made 
to the subdetectors, the trigger system, the run control and the algorithms used in the 
offline reconstruction of tracks, which improved the performance compared to LEP I [11]. 

The major change was the inclusion of the Very Forward Tracker (VFT) [12], which 
extended the coverage of the innermost silicon tracker out to 11° < 6 < 169°. Together 
with improved tracking algorithms and alignment and calibration procedures optimised 
for LEP H, these changes led to an improved track reconstruction efficiency in the forward 
regions of DELPHI. 

A smaller change was the removal of the tungsten nose cone in front of the Small 
Angle Tile Calorimeter, to increase the forward electromagnetic coverage. This change 
had consequences for the luminosity determination discussed in Section 4. 

Changes were made to the electronics of the trigger and timing system which improved 
the stability of the running during data taking [13]. The trigger conditions were optimised 
for LEP H running, to give high efficiency for Standard Model 2- and 4-fermion processes 
and also to give sensitivity to events which might be signatures of new physics. In 
addition, improvements were made to the operation of the detector during the LEP 

^The DELPHI coordinate system is a RH system with the 2 :-axis collinear with the incoming electron beam, the x axis 
pointing to the centre of the LEP accelerator and the y axis vertical. The polar angle 6 is with reference to the 2 :-axis, and 
0 is the azimuthal angle in the y plane, r = \Jx"^ -\- y^. 
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cycle, to prepare the detector for data taking at the very start of stable collisions of 
the e+e“ beams, and to respond to adverse backgrounds from LEP, if they occurred. 
These changes led to an overall improvement in the efficiency for collecting the delivered 
luminosity from ~ 85% at the start of LEP II in 1995 to ~ 95% by the end in 2000. 

During the operation of the DELPHI detector in 2000 one of the 12 sectors of the cen¬ 
tral tracking chamber, the TPC, failed. After September 1®* it was not possible to detect 
tracks left by charged particles inside the broken sector. The data affected corresponds to 
~ 1/4 of that collected in 2000. Nevertheless, the redundancy of the tracking system of 
DELPHI meant that tracks of charged particles passing through that sector could still be 
reconstructed from signals in other tracking detectors. A modihed track reconstruction 
algorithm was used in this sector, which included space points reconstructed in the Barrel 
RICH detector, these points having a signihcant role in the determination of the polar 
angle of tracks. As a result, the track reconstruction efficiency was only slightly reduced 
in the region covered by the broken sector, but on average the resolution on the perigee 
parameters of the tracks was worse than prior to the failure of the sector. The impact 
of the failure of this part of the detector on the different analyses is discussed further in 
Section 5. 


4 Luminosity determination 

The luminosity measurement followed closely the analysis described in [2]. It was 
based on counting the number of Bhabha events in the Small Angle Tile Calorimeter of 
DELPHI, which covered the region between 29 and 185 mrad with respect to the beam 
line. A detailed description of this detector and its performance can be found in [14]. It 
provided a very uniform energy response and an accurate energy resolution of 2.0% for 
100 GeV electrons. The tungsten nose, which was used at LEP I to dehne the inner edge 
of the calorimeter, was removed between data taking in 1995 and 1996. After 1995 the 
dehnition of the geometrical acceptance was entirely based on the reconstructed radii of 
the electron and positron showers. 

The trigger for Bhabha events was prescaled by a factor 3 to 4 to reduce the overall 
trigger rate. A comparison between the measured luminosity and the scalers of the trigger 
shows that the prescaling had no effect on the luminosity measurement [13]. Furthermore 
a prescaled single arm trigger was used to monitor possible trigger inefficiencies, which 
were found to be smaller than 2 x 10“"^. 

In the selection of the Bhabha events, only the most energetic electromagnetic cluster 
in each arm of the STIC was used. To remove the background due to off-momentum 
particles, the following cuts were applied: 

• in each arm, the energy of the cluster was required to be larger than 65% of the 
beam energy; 

• the acoplanarity^ between the two clusters was required to be less than 20°. 

A special trigger, requiring a coincidence between the signal from one arm of the STIC and 
a delayed signal from the other, measured the residual background due to off-momentum 
particles. The measurement showed that it was smaller than 10“^ of the accepted events. 

In order to minimize the sensitivity of the accepted cross-section to the transverse 
position of the interaction point, an asymmetric acceptance was dehned, with a narrow 
side and a wide side. The following cuts were applied to dehne the geometrical acceptance: 

^In general the acoplanarity is defined as the complement of the angle, in the plane transverse to the beam, between 
two tracks. In this case the “tracks” are lines joining the two clusters to the interaction point. 
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Year 

1995 

1996 

1997 

1998 

1999 

2000 

Number of Bhabha events (xlO^) 

156 

144 

384 

1195 

856 

757 


Table 2: Number of Bhabha events at LEP II. 


• the radial position of the reconstructed shower was required to be between 10 and 
25 cm on the narrow side; 

• the radial position of the reconstructed shower was required to be between 9.1 and 
28 cm on the wide side. 

The cuts were chosen at the borders between rings, where the best spatial resolution 
is achieved. The narrow side was alternated between the forward and the backward 
hemispheres in successive events, in order to minimize the sensitivity to the longitudinal 
position of the interaction point. 

The calculation of the accepted cross-section was based on the event generator BH- 
LUMI 4.04 [15], which includes the full 0{a‘^) calculation. The generated events were 
passed through a detailed simulation of the detector, and analysed in the same way as 
the real data. The total accepted cross-section was estimated to be 10.712 ± 0.010 nb, 
at a centre-of-mass energy of 200 GeV. The contribution of the process e’''e“ —77 in 
the selected sample of Bhabha events was calculated to be 0.06%. The theoretical un¬ 
certainty on the estimation of the cross-section is estimated to be ±0.12% [16]. There is 
a factor of two improvement on the error quoted in [3,4]; in the analysis of the results 
and interpretations made in the paper these smaller errors have been applied to all data, 
including those data sets which have not been reanalysed in full. 

The number of selected Bhabhas collected at LEP II energies during the different years 
of data taking are reported in Table 2. 

The experimental systematics in the luminosity measurement were dominated by the 
uncertainty in the radial cut which defines the inner border of the acceptance. The bias 
of the radius reconstruction was carefully studied by comparing the STIC measurement 
with the independent measurement from the two planes of the Silicon Shower Maximum 
Detectors, embedded inside the STIC. The difference between the two measurements was 
monitored throughout data taking and the bias was evaluated to be 250 /rm, correspond¬ 
ing to a systematic uncertainty of ±0.5%. The other sources of uncertainties in the event 
selection were estimated to contribute ±0.04% to the systematics. The luminosity had to 
be corrected for the average displacement of the interaction point from its nominal posi¬ 
tion, as measured by the DELPHI tracking system. The residual systematic uncertainty, 
connected with the £ll-to-£ll £uctuations, was evaluated to be ±0.03%. The uncertainty 
in the beam energy at LEP II propagated into a systematic uncertainty of ±0.04%. A 
detailed list of the contributions to the systematic uncertainty is given in Table 3. The 
overall systematic uncertainty is evaluated to be ±0.5% and it is common to all LEP II 
data taking, except for 1995. For 1995, the presence of the tungsten mask made the 
luminosity determination more accurate, giving an experimental systematic uncertainty 
of ±0.09% on the luminosity determination, for data taken during that year. 
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Source 

AC/C (%) 

Bias on inner edge of the acceptance 

0.5 

Beam energy 

0.04 

Position of the IP 

0.03 

Event selection 

0.04 

Background subtraction 

0.01 

Trigger inefficiency 

0.02 

Total experimental 

0.5 

Total theoretical 

0.12 


Table 3: Contributions to the systematic uncertainty of the luminosity measurement at 
LEP II energies from 1996 onwards. 
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Figure 1: The Standard Model s-channel processes e’''e —ff (upper) and t-channel pro¬ 
cesses in e+e“ —e+e“ (lower). 

5 Analysis of ff final states 

The process e’'“e“ —ff at LEP II is very similar to that at LEP I. In the Standard 
Model the s-channel processes e'''e“ —, e’''e“ —and e’''e“ —qq are described 
by the exchange of virtual photons and Z bosons, shown as Feynman diagrams in Figure 1. 
At LEP I the process is dominated by the formation and decay of the Z. At LEP II, 
the exchange of the photon and the interference between the Z and the photon (which 
is highly suppressed at the Z pole due to the phases of the Z and photon exchange 
amplitudes) are more important. Moreover, the predicted Born level cross-sections for 
the s-channel process are some 2-3 orders of magnitude smaller than at the Z resonance, 
and the forward-backward asymmetries are larger. 

At the peak of the Z resonance QED radiative corrections are signihcant, leading to 
a ~ 30% reduction of the cross-section relative to the Born level prediction. At LEP II 
energies, QED radiative corrections are again signihcant, but here they lead to an increase 
of a factor of 0{2 — 3) in the total cross-section above the Born level predictions. This 
is predominantly through the process of radiative return, in which a photon is radiated 
from the incoming electron or positron, reducing the centre-of-mass energy of the hard 
scattering from the full centre-of-mass energy, ^i, to close to the Z resonance. The 
photon is typically emitted along the direction of the incoming beams, and usually goes 
undetected down the LEP beampipe. The hnal-state fermion pairs in radiative return 
events are, therefore, usually acollinear, though they are typically produced back-to-back 
in the plane transverse to the incoming beams. An important step of the analysis for 
each channel is to determine the reduced centre-of-mass energy of the collisions, . 
This is used to separate events with high collision energies from radiative return events. 
Measured cross-sections and forward-backward asymmetries are quoted for non-radiative 
samples, for which ~ ^/s and for inclusive samples of events, which also include 
radiative return events. The algorithms employed to determine are different in each 
channel. 

The process e’''e“ —e’''e“ is dominated by the t-channel exchange of virtual photons. 
The t-channel processes are shown in Figure 1. At the Z resonance the cross-sections are 
sensitive to s-channel Z exchange. Therefore, at LEP I it was reasonable to subtract off 
the predicted t-channel contributions to measured quantities for e^e“ —e’''e“ collisions in 
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order to extract Z parameters, such as the Z mass and coupling to electrons . Above the Z 
resonance, this subtraction is less appropriate. Physics beyond the Standard Model may 
manifest itself through t-channel processes. No subtraction of the t-channel contributions 
to the measured cross-sections, asymmetries or differential cross-sections has been made. 
A cut on the acollinearity of the electron-positron pair in the hnal state is used to separate 
events with high invariant masses from radiative return events. 

In addition to falling cross-sections for the signal e^e“ —ff processes compared to 
LEP I, the backgrounds from other processes increase at LEP II. The cross-section for 
two-photon collision processes increases as log s and new channels such as e’''e“ —W^W“ 
and e’''e“ —ZZ open at ^/s ~ 2Mw and ^/s ~ 2Mz respectively. 

With some small modihcations to take into account the new radiative return topologies 
and the smaller signal over background ratio, the selection of e’''e“ —ff events at LEP II 
is very similar to that at LEP 1. Classihcation into e’''e“, and qq hnal states is 

based on the the multiplicity of hnal state particles, the responses of the electromagnetic 
and hadronic calorimeters and the muon chambers, and the momenta of charged particles 
measured in the tracking system. 

Event simulation 

To determine selection efficiencies and backgrounds for each analysis, events were 
simulated using a variety of event generators and the DELPHI detector simulation [11] 
and passed through the full data analysis chain. To allow studies of the data taken after 
September 1®* 2000, samples of events were simulated dropping information from the 
broken sector of the TPC. 

The s-channel e’''e“ —ff processes were simulated with KK 4.14 [17], while events in 
the e’''e“ —e+e” channel were simulated with BHWIDE 1.01 [18]. The fragmentation of 
qq events into hadrons was simulated using PYTHIA 6.156 [19,20], with parameters tuned 
to DELPHI data [21]. Spin effects in r-lepton decays were handled by TAUOLA 2.6 [22] 
using the helicity approximation. There were typically 1 million simulated events per 
energy for the e’'“e“ —> e’'“e“ and e’''e“ —qq channels, with 100,000 simulated events for 
the e'''e“ —> channel and 200,000 simulated events for the e’''e“ —channel at 

each energy. 

Four-fermion background events, including high-mass two-photon collision events, were 
simulated with the generator WPHACT 2.0 [23]. The generated events were divided into 
3 samples [24]; the hrst dominated by charged current processes, e+e“ —W"*"W“; the sec¬ 
ond dominated by neutral current processes with topologies similar to e’''e“ —ZZ events; 
and the third sample being neutral current samples dominated by multi-peripheral pro¬ 
cesses, e'''e“ —> e’''e“ff. Low mass two-photon collisions were simulated with BDKRC [25] 
for leptonic hnal states and with PYTHIA for hadronic hnal states. Since the samples of 
4-fermion events were generated with certain kinematic cuts, the background subtraction 
involved an extrapolation to estimate the backgrounds coming from events which were 
not simulated. Theoretical uncertainties amount to ±5% on the e’''e“ —e’''e“ff samples, 
±0.5% on the e’''e“ —W'''W“ and ±2% on the e’''e“ —ZZ samples. Where possible 
the real data were used to control the simulation predictions for the backgrounds in the 
selected samples. 

Theoretical predictions and signal definition 

The measurements reported in this paper are compared to theoretical predictions, from 
the Standard Model and from models which include physics beyond the Standard Model. 
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Throughout this paper Standard Model predictions are taken from the semi-analytical 
QED corrected, Improved Born Approximation computations of ZFITTER [26] for 
e'''e“ —1 e+e“ —r+r“ and e’''e“ —qq, and from BHWIDE [18] for e+e“ —e+e“. 
ZFITTER version 6.36 ^ was used with the following central values for input parameters 


Mz 

= 91.1875 GeV/c^ 

mt 

= 175.0 GeV/c^, 

Mh 

= 150.0 GeV/c^ 

<As(Mz) 

= 0.118, 


= 0.02761, 


where is the contribution to the running of the electromagnetic coupling constant, 

a, due to contributions from hadronic loops containing 5 quark flavours. To make the 
comparison it was necessary to match the signal dehnitions in the data, simulation and 
semi-analytical computations. There were several choices to be made: dehnition of ; 
subtraction or inclusion of QED corrections from the interference between Initial State 
Radiation (ISR) and Final Sate Radiation (FSR); inclusion of ff events with additional 
radiated pairs of fermions; angular acceptance and acceptance. The signal dehnition 
adopted here was chosen to make analysis of the non-radiative events as straightforward 
as possible: 

• for e^e“ —and e’''e“ —> is taken to be the invariant mass of the 

fermion pair, Mfp For e’''e“ —qq, is taken to be the invariant mass of the 
fermion pair with any FSR included - for s-channel processes this corresponds to 
computing the invariant mass, Q, of the virtual propagator; 

• QED corrections from the interference between ISR and FSR are included for 
e’''e“ —and e’''e“ —. These corrections are included in the simulated 
events. For e’''e“ —qq these corrections are not included in the simulated events, 
furthermore, the dehnition = Q is ambiguous in the presence of interference 
between ISR and FSR. Therefore, for the e’''e“ —qq channel, the signal is dehned 
as having no ISR-FSR interference. Corrections are applied to the data to remove 
the ehects of the interference; 

• events with additional radiated fermion pairs are subtracted as part of the 4-fermion 
background. The bulk of the pairs come from internal conversion of ISR photons into 
e'''e“ pairs, which are lost in the beampipe, and which are topologically equivalent 
to radiative return events. However, for the non-radiative samples this background 
is small compared to the signal; 

• the total cross-sections and forward-backward asymmetries are quoted in the full dvr 
acceptance for e’''e“ —, e+e“ —r+r“and e’''e“ —qq. For e’''e“ —and 
e’''e“ —> T^T~ this involves an extrapolation from the hducial volume of the detector 
using generated events. For the electrons the measurements are made within the 
acceptance 44° < 6 < 136° - which corresponds to the acceptance of the HPC. 
Measurements in the FEMC region, 12° < 6 < 35° and 145° < 6 < 168°, are not 
reported here, since they are subject to greater experimental uncertainties and the 

=*The following ZFITTER flags were used: AFBC: 1 SCAL: 0 SORE: 0 AMT4: 4 BORN: 0 BOND: 2 CONV: 2 FOT2: 3 
GAMS; 1 DIAG; 1 BARB; 2 PART: 0 POWR: 1 PRNT; 0 ALEM: 2 QGDC: 3 VPOL: 1 WEAK; 1 FTJR: 1 EXPR; 0 
EXPF: 0 HIGS: 0 AFMT; 1 CZAK: 0 PREC:10 HIG2: 0 ALE2: 3 GFER: 2 ISPP: 2 FSRS: 1 MISC: 0 MISD: 1 IPFC: 5 
IPSG: 0 IPTO:-l FBHO: 0 FSPP: 0 FUNA: 0 ASCR: 1 SFSR: 1 ENUE: 1 TUPV: 1. For e+e" ^ qq FINR; 0 INTF: 0, 
while for e+e“ —> and e+e“ —> FINR: 1 INTF: 2 For non-radiative samples FSRS: 0, while for inclusive 

samples FSRS; 1. 
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cross-sections are dominated by the t-channel process. The differential cross-sections 
for e+e“ —and e’''e“ —are quoted within the hducial volume of the 
detector, with a cut at | cos6*| = 0.97 for e’''e“ —and 0.96 for r’''r“final states; 

• for the inclusive samples the acceptances for t~^t~ hnal states are 

\/i'>75 GeV and V^>0.10 ^/s for qq final states. For the non-radiative samples 
a/s'/v^> 0.85 for all these processes. For e’''e“ hnal states the cut on is replaced 
by a cut on the acollinearity, 9acoi < 20°; both s and t channel processes leading to 
e’''e“ hnal states are considered as signal. 

In the following sections the analyses of the diherent hnal states are discussed and 
results of the measurements are presented. In all cases the latest theoretical predictions 
have been used, updating values given in previous publications [3,4]. 

Experimental uncertainties 

Although a large number of possible biases and sources of uncertainty were investigated 
for all hnal states, only those sources of bias and uncertainty which lead to signihcant 
systematic errors for each particular analysis are described in the sections below. Sources 
of bias and uncertainty which are negligible for a particular channel are not described. 

For measurements of the diherential cross-sections for the e’''e“ —and 
e’''e“ — T~^T~ channels, some of the bins used for the analysis contain only a small number 
of observed events. For these measurements statistical errors were computed both from 
the square root of the number of observed events and of the number of events expected 
from the Standard Model. The second error provides a suitable weight which can be used 
to combine measurements from diherent energies and can also be used when combining 
data from diherent LEP experiments or htting small deviations from the Standard Model 
to the measurements. In all other cases the statistical errors were computed solely from 
the square root of the number of observed events. 

Analyses of the various final states 

5.1 final states 

An analysis of e+e“ hnal states at e’''e“ collision energies of 189 GeV and above 
is presented. This updates the analysis of data taken at 189 GeV [4]. Gompared to 
the previous analysis, an explicit correction is applied for charge misidentihcation in 
the measurement of the forward-backward asymmetry. New results are presented for 
diherential cross-sections. Data taken at collision energies below 189 GeV [3,4], have not 
been reanalysed. 

5.1.1 Analysis 
Run Selection 

Runs were selected requiring that the Vertex Detector (VD), the TPG and the elec¬ 
tromagnetic calorimeters in the barrel of DELPHI (HPG) were operative. For the 2000 
data fraction with a TPG sector dead, the requirement for the TPG was restricted to the 
remaining sectors of the detector. 
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Event selection 

Events were selected with two almost-independent sets of experimental cuts, chosen 
in such a way as to minimize the correlations between the two sets [1,2]. Only the barrel 
region of DELPHI was used for this analysis. In each selection, both the electron and the 
positron were required to be within the range 44° < 6 < 136° and the acollinearity was 
required to be smaller than 20°. A cut in polar angle at 90 ± 2° was applied to remove 
the region with neither TPC nor HPC coverage. 

The hrst set of cuts (selection A) is based on requirements of calorimetric-energy 
clusters associated to hits in the Vertex Detector. In particular events were selected if 
they had: 

• at least two electromagnetic clusters in the HPC one with energy above 75% of the 
beam energy and another above 55% and with an acollinearity angle between the 
clusters less than 20°; 

• at least two track segments in opposite hemispheres seen by the VD and no more 
than four in total; events with a 2-versus-2 topology were excluded; 

• no energy in the last two layers of the hadron calorimeter (HCAL) in correspondence 
with electromagnetic clusters detected at large distance (> ±1° in 0) from the HPC 
sector boundaries. 

The second selection (selection B) is based on the charged-particle momentum as 
measured by the tracking system. In particular events were selected if they had: 

• at least 2 charged-particle tracks, of momentum greater than 1.5 GeV/c and distance 
of closest approach to the nominal vertex position less than 5 cm, seen by the 
DELPHI tracking system outside the VD, with acollinearity less than 20° and no 
more than four tracks in total; the 2-versus-2 track topology was excluded; 

• the quadratic sum {prad) of the two momenta of the highest-momentum charged 
particle in each hemisphere greater than 0.99y^/2; 

• no energy observed in the last three layers of HCAL associated to the impact points 
of the two highest-momentum charged particles; 

• the OD hit pattern associated to the impact points of the tracks compatible with 
the pattern of a particle showering in or before the OD, or giving back-scattering 
from the HPC calorimeter; 

• no tracks identified as a muon. 

Estimation of the selection efficiency and background 

Considering the selections A and B as independent, the efficiency of each of them 
could be easily computed by a comparison of the number of events selected by each one 
separately or by both simultaneously. This is much easier at LEP H energies, given the 
t-channel dominance, compared to the analysis at the Z resonance because of the much 
smaller background ( r events estimated from simulation) as shown in Figure 2(left). 

The efficiency of the two selections at the different energy points is given in Table 4. 
In both selections the measured efficiencies did not include the loss due to the exclusion 
of the polar angle region around 90°. The total number of events selected by selection A 
at the different centre-of-mass energies is shown in Table 4. 

The simulated e’''e“ events were used to estimate and remove the bias caused by the 
correlation between the two selections due to the detector structure or to the kinematics 
of the events. The bias on the efficiency was found to be (0.4 ± 0.7)%. However in 
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Figure 2: The energy of the two most energetic clusters (left) and the uncorrected cos 9 
distribution (right). Both distributions correspond to data taken at ^/s ^-^200 GeV and 
the left plot does not include the energy cuts (indicated by the arrows) applied in selection 
A. The histogram corresponds to the simulated events and the dark area shows the r 
background. 








































14 




e'''e 

—>■ e'''e 






Energy point (GeV) 

189 

192 

196 

200 

202 

205 

207 

Energy (GeV) 

188.63 

192.17 

196.10 

200.12 

202.07 

204.88 

206.59 

Lumi (pb“^) 

155.11 

25.12 

76.16 

83.07 

40.05 

88.55 

128.39 

No. Events (A) 

3179 

518 

1568 

1554 

778 

1500 

2126 

Efficiency (A) (%) 

91.9 

92.3 

93.7 

92.8 

93.6 

92.6 

92.6 

Efficiency (B) (%) 

77.4 

76.6 

76.4 

76.7 

75.1 

72.6 

72.6 

Background (%) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 


Table 4: Details of LEP II analysis for e'^e~ e’''e“ channel. The table shows the actual 

centre-of-mass energy and luminosity analysed at each energy point, the number of events 
selected by selection A and the efficiencies of selecting events with selections 

A and B and the background selected with selection A. 

the 2000 data new algorithms were used for track reconstruction based on VD hits. An 
increased correlation between the two selections was expected. In order to measure it, 

1999 data reconstructed with the new and the old algorithm were compared. This gave a 
correction factor of 1.016 ± 0.003 to be applied to the cross-section measurements of the 

2000 data. 

5.1.2 Results 
Cross-sections 

The cross-section was obtained as: 

O-ee = 7- Cf + (1) 

LCa 

where Na is the number of events selected by selection A and ca is the efficiency of this 
selection, C is the luminosity, c/ is a correction factor including the bias in the evaluation 
of €a due to the correlation between the two selection methods and also the background 
subtraction, and <7®° is the correction for the central region (90±2)°. The value of <7®° was 
computed by TOPAZO [271 and ALIBABA [281 and its value ranges from 0.43 ± 0.04 pb 
to 0.36 ±0.04 pb. 

The statistical error on the cross-section includes the statistical uncertainty on the de¬ 
termination of the efficiency, taking into account the statistical uncertainties on correlated 
and uncorrelated subsamples from selections A and B. 

All the values of the measured cross-sections are given in Table 5 and shown in Fig¬ 
ure 10. 

Forward-backward asymmetries 

In order to measure the forward-backward asymmetry the charge of the event was 
dehned as positive when the positron was in the forward hemisphere with respect to 
the incoming positron direction, and negative in the opposite case. In the e’''e“ —e’''e“ 
events, in addition to the canonical charge dehnition from reconstructed tracks, it is 
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e'*’e —^ e'*’e 

Vi 

(GeV) 

Oacol < 20° 

Vi 

(GeV) 

Oacol < 20° 

O-ee (pb) 

A® 

^FB 

O-ee (pb) 

A® 

^FB 

130 

42.00 ±4.00 ±0.78 

(48.70) 

0.810 ± 0.060 ± 0.003 

(0.810) 

192 

22.71 ± 0.98 ± 0.24 

(22.13) 

0.831 ± 0.024 ± 0.003 

(0.820) 

136 

47.10 ±4.20 ±0.73 

(44.60) 

0.890 ± 0.040 ± 0.003 

(0.810) 

196 

22.33 ±0.55 ±0.23 
(21.24) 

0.823 ± 0.014 ± 0.003 
(0.821) 

161 

27.10 ± 1.80 ±0.43 

(31.90) 

0.820 ± 0.040 ± 0.003 

(0.830) 

200 

20.52 ±0.51 ±0.21 

(20.36) 

0.788 ±0.016 ±0.003 

(0.823) 

172 

30.30 ± 1.90 ±0.45 

(28.00) 

0.810 ± 0.040 ± 0.003 

(0.830) 

202 

21.11 ± 0.74 ± 0.22 

(19.97) 

0.831 ± 0.020 ± 0.003 
(0.822) 

183 

25.63 ±0.76 ±0.26 
(24.54) 

0.814 ±0.017 ±0.003 

(0.817) 

205 

18.94 ± 0.48 ± 0.21 

(19.33) 

0.797 ±0.016 ±0.004 

(0.820) 

189 

22.73 ±0.40 ±0.23 

(22.93) 

0.804 ± 0.010 ± 0.003 

(0.820) 

207 

18.52 ±0.39 ±0.20 

(19.07) 

0.820 ±0.012 ±0.004 
(0.822) 


Table 5: Measured eross-sections and forward-baekward asymmetries for non-radiative 
events. The statistieal error is followed by the total systematie error. In 
parentheses the expected values as computed by BHWIDE, which have a precision of 
±2% on (Tee ±0.02 on are given. Results are quoted for an acceptance of 

44° < 9 <130°. 
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possible to correlate the charges of an event by looking at the effects of the bending dne 
to the magnetic field on the impact position of HPC clusters, giving a high redundancy 
on the charge determination. The latter method was used to determine the charge of an 
event in cases where the reconstructed charges of the tracks were equal. 

Given the high expected asymmetry the measurements have to be corrected for the 
residual wrong charge assignments. The asymmetry was corrected by a factor l+2,^± with 
= (7.2 ± 1.4) ■ 10“^ corresponding to the charge misassignment which was determined 
from simulated events: 


4° - 

^FB — 


N+-N_ 
N+ + N_ 


■(l + 2e±), 


( 2 ) 


where iV+ and iV_ are the number of events with positive and negative charge, respec¬ 
tively. In the previous publications this correction was not applied and an asymmetric 
error was given. The forward-backward asymmetry is corrected also for the missing 
central region; 
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( 3 ) 


where and are the computed cross-sections for the regions [88°, 90°] and [90°, 92°] 
respectively. 

All the values of the measured forward-backward asymmetries are given in Table 5 
and shown in Figure 11. 


Differential cross-sections 


As with the measurement of the forward-backward asymmetry, measurements of the 
differential cross-sections require the determination of the charge of the electron and 
positron, and a correction has to be applied for mismeasurements of the charge. A bin 
by bin correction was applied to the events detected in each polar angle bin which 
was largest in the very forward and backward bins due to the large asymmetry in the 
cross-sections: 

/VJ = /V„.-(l+{±)(4) 

where Ng. is the number of events detected in the opposite hemisphere. The bin —0.72 < 
cos 6 < —0.54 has the largest correction which is about 25 to 30%. The differential 
cross-section is obtained as 


d(j/d cos 6*^ 




Ceg. ■ A cosO' 


( 5 ) 


where ee^ is a bin correction factor dominated by the efficiency evaluation and A cos 9 
is the bin width. The 90° region correction is taken into account by adding / A cos 9 
and a'^/A cos 9 to the corresponding bins. The uncorrected cos 9 distribution is shown 
in Figure 2(right). The differential cross-sections for the 1998, 1999 and 2000 data are 
given in Table 6. 


Systematic errors 

Apart from the luminosity, systematic uncertainties arise from the event selection 
correlation, acceptance dehnition and from the background subtraction. The largest un¬ 
certainty comes from the selection correlation (±0.7%) due to the statistics of simulated 
events. The simulated data for 2000 were produced exclusively with a new tracking 
algorithm giving large correlations (which were measured directly from the data) and 
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P"*" P 


yCg ~ 196 


dcr/dcos^ (pb) 

COS 9 

SM 

Measurement 

[-0.72,-0.54] 

1.68 

1.73±0.41±0.09 

[-0.54,-0.36] 

2.21 

2.37±0.43±0.04 

[-0.36,-0.18] 

2.77 

3.03±0.48±0.03 

[-0.18, 0.00] 

3.89 

4.13±0.50±0.23 

[ 0.00, 0.09] 

5.52 

7.23±0.75±0.34 

] 0.09, 0.18] 

6.97 

7.01±1.04±0.07 

] 0.18, 0.27] 

9.53 

11.27±1.23±0.09 

] 0.27, 0.36] 

13.38 

14.29±1.47±0.13 

] 0.36, 0.45] 

19.39 

23.39±1.72±0.19 

] 0.45, 0.54] 

29.77 

30.14±2.14±0.29 

] 0.54, 0.63] 

48.46 

51.75±2.78±0.48 

] 0.63, 0.72] 

81.88 

81.55±3.51±0.81 


yS ~ 192 


dcr/dcos^ (pb) 

COS 6 

SM 

Measurement 

[-0.72,-0.54] 

1.81 

2.03±0.73±0.09 

[-0.54,-0.36] 

2.29 

1.87±0.80±0.04 

[-0.36,-0.18] 

2.93 

2.40±0.87±0.03 

[-0.18, 0.00] 

4.05 

5.36±1.00±0.23 

] 0.00, 0.09] 

5.64 

4.05±1.40±0.34 

] 0.09, 0.18] 

7.44 

11.39±2.18±0.08 

] 0.18, 0.27] 

9.83 

9.74±2.20±0.10 

] 0.27, 0.36] 

13.80 

12.03±2.64±0.14 

] 0.36, 0.45] 

20.35 

20.72±3.34±0.21 

] 0.45, 0.54] 

31.32 

29.61±3.88±0.32 

] 0.54, 0.63] 

50.21 

54.58±4.79±0.51 

] 0.63, 0.72] 

85.17 

87.93±6.28±0.87 


yCs ~ 189 


dcr/dcos^ (pb) 

cos 6 

SM 

Measurement 

[-0.72,-0.54] 

1.89 

2.69±0.31±0.09 

[-0.54,-0.36] 

2.36 

2.31±0.32±0.04 

[-0.36,-0.18] 

2.98 

2.53±0.34±0.03 

[-0.18, 0.00] 

4.23 

4.73±0.36±0.23 

] 0.00, 0.09] 

5.86 

5.27±0.49±0.34 

] 0.09, 0.18] 

7.55 

8.18±0.77±0.07 

] 0.18, 0.27] 

10.16 

9.79±0.88±0.10 

] 0.27, 0.36] 

14.44 

12.58±1.06±0.14 

] 0.36, 0.45] 

21.15 

21.48±1.31±0.21 

] 0.45, 0.54] 

32.75 

31.21±1.60±0.32 

] 0.54, 0.63] 

52.28 

52.76±2.04±0.51 

] 0.63, 0.72] 

87.65 

87.85±2.59±0.86 


y/s ^ 200 


dcr/dcos^ (pb) 

cos 9 

SM 

Measurement 

[-0.72,-0.54] 

1.57 

1.73±0.38±0.08 

[-0.54,-0.36] 

2.07 

2.81±0.40±0.03 

[-0.36,-0.18] 

2.66 

2.86±0.45±0.03 

[-0.18, 0.00] 

3.73 

4.95±0.48±0.23 

] 0.00, 0.09] 

5.23 

5.23±0.73±0.34 

] 0.09, 0.18] 

6.72 

6.26±0.95±0.07 

] 0.18, 0.27] 

9.02 

9.14±1.17±0.09 

] 0.27, 0.36] 

12.85 

13.16±1.37±0.13 

] 0.36, 0.45] 

18.57 

20.06±1.64±0.18 

] 0.45, 0.54] 

29.03 

26.43±2.03±0.29 

] 0.54, 0.63] 

46.21 

46.00±2.58±0.46 

] 0.63, 0.72] 

78.54 

77.96±3.34±0.78 


205 


dcr/dcos^ (pb) 

COS 9 

SM 

Measurement 

[-0.72,-0.54] 

1.58 

1.94±0.37±0.07 

[-0.54,-0.36] 

1.96 

2.21±0.38±0.03 

[-0.36,-0.18] 

2.56 

3.01±0.43±0.03 

[-0.18, 0.00] 

3.57 

3.73±0.44±0.23 

] 0.00, 0.09] 

4.88 

5.21±0.59±0.34 

] 0.09, 0.18] 

6.33 

5.48±0.98±0.07 

] 0.18, 0.27] 

8.49 

8.49±1.08±0.09 

] 0.27, 0.36] 

12.27 

11.53±1.30±0.13 

] 0.36, 0.45] 

17.80 

17.62±1.56±0.19 

] 0.45, 0.54] 

27.45 

28.17±1.94±0.29 

] 0.54, 0.63] 

44.01 

42.14±2.45±0.46 

] 0.63, 0.72] 

74.19 

70.86±3.14±0.77 


y/s ^ 202 


dcr/dcos^ (pb) 

cos 9 

SM 

Measurement 

[-0.72,-0.54] 

1.61 

1.50±0.55±0.09 

[-0.54,-0.36] 

2.03 

1.96±0.57±0.04 

[-0.36,-0.18] 

2.58 

2.98±0.62±0.03 

[-0.18, 0.00] 

3.68 

3.40±0.64±0.23 

] 0.00, 0.09] 

4.99 

4.55±0.86±0.34 

] 0.09, 0.18] 

6.63 

6.16±1.39±0.07 

] 0.18, 0.27] 

8.90 

5.82±1.57±0.09 

] 0.27, 0.36] 

12.44 

13.02±1.93±0.13 

] 0.36, 0.45] 

18.44 

18.85±2.40±0.19 

] 0.45, 0.54] 

28.11 

33.61±2.85±0.28 

] 0.54, 0.63] 

45.70 

48.97±3.74±0.46 

] 0.63, 0.72] 

76.92 

84.69±4.74±0.78 


yi ~ 207 


dcr/dcos^ (pb) 

COS 9 

SM 

Measurement 

[-0.72,-0.54] 

1.53 

1.25±0.30±0.07 

[-0.54,-0.36] 

1.95 

2.25±0.32±0.03 

[-0.36,-0.18] 

2.55 

2.11±0.35±0.03 

[-0.18, 0.00] 

3.41 

3.85±0.36±0.23 

] 0.00, 0.09] 

4.84 

5.84±0.49±0.34 

] 0.09, 0.18] 

6.34 

6.30±0.81±0.07 

] 0.18, 0.27] 

8.30 

7.59±0.89±0.09 

] 0.27, 0.36] 

11.82 

10.86±1.06±0.12 

] 0.36, 0.45] 

17.70 

16.84±1.29±0.19 

] 0.45, 0.54] 

26.95 

25.22±1.60±0.28 

] 0.54, 0.63] 

43.65 

42.40±2.03±0.46 

] 0.63, 0.72] 

73.37 

72.69±2.60±0.77 


Table 6: Differential cross-sections for non-radiative e’''e“ events at cen- 

tre-of-mass energies from r^l89 to 207 GeV. The tables show the bins and the predicted 
Standard Model differential cross-sections and the measurements with statistical and ex¬ 
perimental systematic errors. 
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e’^e — 

e'^e 


Source 

Aa/a 

(non-rad.) 

AApB 

(non-rad.) 

Efficiency 

76 

- 

Acceptance 

48 

20 

Charge misidentification 

- 

20 

Backgrounds 

20 

20 

Luminosity 

58 

- 

TPC sector instability 

- 

20 

Total uncertainty 

109 

35 

Correlated 

102 

35 

Uncorrelated 

38 

- 


Table 7: Systematic uncertainties on the measured cross-sections and forward-backward 
asymmetries for data taken at y/s ^201 GeV. All numbers in units of 10~^. The total 
uncertainty does not include the error due to TPC sector instability which applies only 
to a part of the 2000 data. The correlated error component includes errors correlated 
between energies and channels and those correlated with other LEP experiments. 
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therefore the statistical error cannot be reduced. The error on the cross-section due to 
the acceptance dehnition arises from the uncertainty on the absolute polar angle determi¬ 
nation. For events having fully reconstructed charged particle tracks this uncertainty is 
very small (< 0.02°), while for those having the acceptance dehned by HPC clusters it is 
larger (about 0.04°). The latter sample has an additional contribution due to the poorer 
precision (±0.2°) on the cluster polar angle determination. The total contribution to the 
acceptance definition depends on the fraction of events with HPC-based acceptance and 
it ranges from ±0.24% to ±0.30%. Beam energy spread effects on the acceptance were 
investigated and they were found to be ±0.28%. 

The corrections applied for the ±2° polar angle hducial cut around 90° in the anal¬ 
ysis of the total and differential cross-sections were computed at the different energies 
by using the program TOPAZO and checked with ALIBABA. No signihcant differences 
were found between the two generators. An error of ±0.25% was assigned to the total 
cross-section correlated between all energies. The rate of events with two same-charged 
particle tracks was compatible between data and simulation and consequently the charge 
misassignment was determined from simulated events. The charge misassignment error 
was determined by the statistics of the simulated events and it amounts to ±0.002 on the 
forward-backward asymmetry, correlated among all energies. 

The TPC failure in the last part of the 2000 run had a clear influence on the Selection B 
efficiency: 63.0% instead of 75.5% obtained during the first part of the 2000 run. However 
the overall effect on the cross-sections and asymmetries was not very signihcant: the sta¬ 
tistical error increase was of the order of a few percent and only an additional systematic 
error of ±0.3% on the asymmetry was estimated because of the charge determination. 

A breakdown of the systematic errors for data taken at ^/s ~ 207 GeV is given in 
Table 7. 

5.2 final states 

An analysis of hnal states at e’''e“ collision energies of 183 GeV and above is 

presented. Data taken at energies of 183 GeV and 189 GeV [4] have been reanalysed to be 
consistent with data taken at higher energies. Gompared to the previous analysis, the new 
analysis benehted from an increase in Monte Garlo statistics available in the whole LEP 
II data set. This made it possible to perform detailed comparisons of data and simulated 
events, which ultimately led to improvements in the estimation of the systematic errors on 
the measurements. For data taken at collision energies below 183 GeV [3], the statistical 
errors are so large that improvements in the systematic errors would be negligible. These 
data have not been reanalysed. 

5.2.1 Analysis 
Run selection 

Runs in which important components of the detector, such as the muon chambers, 
were not functioning satisfactorily were removed from the analysis. 

Event selection 

Events in the channel e’''e“ —were selected by first imposing the following 
kinematic requirements: 

• at least two reconstructed tracks from charged particles; 
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• to allow for photon conversions or splitting of tracks, events with up to 7 recon¬ 
structed tracks were accepted; 

• the measured momenta of the two highest-momentum reconstructed charged parti¬ 
cles had to be greater than 15 GeV/c; 

• the polar angles of the two highest-momentum reconstructed charged particles had 
to lie between 14° and 166°; 

• to suppress backgrounds from e’''e“ —events, the momentum of the third 
highest-momentum charged particle had to be less than 5 GeV/c in all cases where 
the energy of the leading charged particle was less than 90% of the beam energy. 

Individual charged particles were identified as muons if: 

• there was at least one hit in the muon chambers associated to the reconstructed 
track; or 

• the associated energy deposits in the HGAL were consistent with a minimum ionising 
particle; or 

• the associated energy in the electromagnetic calorimeters was less than 1.5 GeV. 

The two highest-momentum charged particles had to be identihed as muons. Events were 
rejected if; 

• either of the two highest-momentum charged particles had no associated muon cham¬ 
bers hits and the energy per layer of the HGAL exceeded 5 GeV; or 

• the electromagnetic energy associated to either of the two highest-momentum par¬ 
ticles exceeded 10 GeV and neither track had associated muon chamber hits - this 
removed e^e“ —> e^e“ events. 

To reduce cosmic ray backgrounds, the reconstructed tracks of the two highest-momentum 
charged particles had to appear to come from the interaction region: 

• for tracks which had associated VD hits, the modulus of the impact parameter 
relative to the beamspot in the plane transverse to the beam axis had to be less 
than 0.1 cm, for at least one of the tracks; 

• for tracks without hits in the VD, the modulus of the impact parameter had to be 
less than 1.0 cm, for both tracks, reflecting the worse resolution for these tracks; 

• in all cases, the impact parameter with respect to the beamspot in the direction 
parallel to the beam axis had to be less than 2.0 cm, this cut being set by the length 
of the electron and positron bunches. 

\/~P reconstruction 

Having selected suitable e’''e“ —events, kinematic hts were used to determine 
the invariant mass of the pair, The kinematic hts attempted to match the 

measured event to several possible event topologies, using energy and momentum con¬ 
straints to improve the determination of the 3-momenta of the reconstructed muons. Use 
of the kinematic hts improved the resolution of the measured invariant masses of the 
pairs, compared to using the directly measured momenta of all muons, or using the 
measured angles of the muons, assuming a single photon emitted down the beampipe. As 
well as using the tracks selected as muons, up to one reconstructed photon per event was 
included in the ht. Electromagnetic energy clusters, not associated to any reconstructed 
track, were considered as photons if they had an energy exceeding 5 GeV. 

The following kinematic hts were tried in the order given: 
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DELPHI 



Figure 3: The distribution of d^i + 6^2 — 18(F for the two highest momentum muons in 
data and simulation. For back-to-back events the sum should be (T. This is broadened by 
radiation, resolution and beam energy spread. Data from 1997-2000 are shown. 

• if a photon was found, the event was htted assuming a seen photon and a possible 
additional unseen photon; 

• a £t in which there was assumed to be no photons radiated in the event; 

• a fit assuming a single photon emitted along the beampipe; 

• a fit assuming a single photon emitted in any direction but unseen in the detector. 

In each case, if the £t gave a satisfactory the invariant mass of the pFyi~ pair was cal¬ 
culated from the fitted momenta of the particles. Otherwise, the next £t was attempted. 
If none of the fits gave an acceptable the measured momenta of the muons were used 
to determine the invariant mass of the p~ pair. Events were divided into inclusive and 
non-radiative classes according to the dehnitions given in the introduction to Section 5. 

In the non-radiative class, migrations into and out of the sample were determined from 
the simulated events. The migrations amount to a correction of between 1.6% and 3.6% 
on the measured cross-sections. 

Estimation of the selection efficiency and background 

Selection efficiencies and backgrounds were determined from simulated events, but 
critical components of the analysis, such as the efficiency for reconstructing and identify¬ 
ing muons, were determined from the data. To obtain good agreement between the data 
and the predictions of simulations it was necessary to apply a number of corrections. 

Using e^e“ —>• p'^p~ events, collected at a/s ~ Mz at various times during 1997 
through to 2000, smearings and shifts in the mean values for distributions of 1 /p were 
computed for different charges of particle, different ranges of cosO, and different com¬ 
binations of tracking detectors involved in the track reconstruction. These corrections 
were applied to data taken at LEP II energies, improving the agreement between data 
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and simulations in the distribution of the momenta of particles at these higher centre-of- 
mass energies. Application of these corrections led to only small changes in the expected 
numbers of events selected in the inclusive and non-radiative classes. This reflects the 
fact that the relative weights given to measurements of angles in the kinematic hts were 
larger than the weights given to the less precisely measured momenta of particles. 

In the simulations of e+e” —events, the momenta of the incoming electrons 
and positrons were taken to be equal and opposite. However, in the data the particles 
in the incoming bunches of electrons and positrons had a momentum spread of ~ 0.2 
GeV/c, such that the energies of the electron and positron before collision were not 
necessarily equal. This led to an observable broadening of the acollinearity distribution 
at low acollinearities, in the data. By applying boosts to the simulated events, 

after the full detector simulation, reflecting the momentum spread of the real beams, 
this effect was taken into account. This led to improved agreement between data and 
simulation in distributions which were used to study the performance of the kinematic 
hts. The agreement between data and simulation after inclusion of the appropriate boosts 
is shown in Figure 3. Due to the relatively loose cuts in v^/v^, the application of the 
boosts in the simulation made a relatively small change to the expected number of events 
selected in the non-radiative sample. 

Corrections were also applied to bring agreement between the data and simulation 
for the numbers of muon chamber hits associated to reconstructed tracks. These were 
determined from the high energy data themselves. Part of the corrections were to take 
into account unsimulated high voltage trips of the chambers which actually happened 
during the collection of the LEP II data set. Because the muon identihcation criteria 
select events from a union of 3 highly overlapping samples, application of the correction 
had only a limited impact on the expected numbers of particles identihed as muons. 

After applying the corrections above, differences between data and simulation were 
still found in the efficiency to reconstruct and identify single muons. High purity samples 
of single muons with high momentum were selected, and the efficiency to reconstruct and 
identify a single muon was determined from the number of times a second muon was found 
accompanying the hrst muon. The backgrounds from events other than e+e“ — 
events, such as e+e“ —r+r“ and e’''e“ —which could lead to only one muon 
and one other track in the detector, were negligible. Averaging over all data from ^/s ~ 
183 GeV and above, it was found that the efficiency for selecting e’''e“ —events 
determined from the simulation had to be reduced by (1.08 ± 0.49)%. The error on this 
value, which is dominated by statistics of events selected in data, is the most signihcant 
systematic error for the measurement of the non-radiative cross-sections, and is correlated 
between measurements at all energies. 

The efficiency for selecting e’''e“ —events determined from the simulation was 
also corrected for the efficiency for triggering DELPHI on these events. This was esti¬ 
mated using the redundancy of the DELPHI trigger system in selected e’''e“ — 
events. It was found to be stable over all years of data taking. Averaging over all years 
the efficiency was determined to be (99.82 ±0.07)%. No signihcant variation in efficiency 
was found over cos 6 within the acceptance of the analysis. 

Backgrounds from e’'“e“ —e’''e“ —> W’''W“, e’'“e“ —ZZ and e’'“e“ —> e’''e“ff 
processes were estimated from samples of simulated events. Backgrounds from 
e'''e“ —> e'''e“, e+e” —qq and two-photon collisions were estimated to be negligible. The 
background from cosmic rays inside the selected sample was estimated by extrapolating 
the numbers of cosmic rays failing the cosmic ray rejection cuts, from within a looser selec- 



23 


e+e — p.'^p, 

Background 

Non-radiative 

Inclusive 

source 

background 

background 

e'''e“ —> r+r“ 

0.08 ± 0.02 

0.94 ± 0.03 

e+e- ^ W+W- 

0.16 ± 0.07 

3.56 ± 0.21 

e+e- ^ ZZ 

0.26 ± 0.06 

1.01 ± 0.06 

e+e“ —»• e+e“ff 

0.10 ± 0.06 

3.75 ± 0.59 

Cosmic rays 

0.44 ± 0.04 

0.45 ± 0.04 

Total 

1.04 ± 0.12 

9.72 ± 0.63 


Table 8: Residual background levels for data taken at ^/s ^207 GeV (in % relative to the 
number of selected e+e“ —>• candidates). 
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Figure 4: The reconstructed y/R/y/s distribution (left) and the reconstructed yfR distri¬ 
bution (right) for e’''e“ — events. Both plots show data from 1997-2000. In the 
right hand plot all selection criteria except for the cut on the pGinvariant mass were 
applied. 

tion on the impact parameters of the reconstructed tracks into the selected e''“e“ —> pGpi~ 
region. The estimated backgrounds are given in Table 8. 

The distributions of /y/s for events in the non-radiative region and y/^ in the 
radiative return region are shown in Figure 4. There is good agreement between data 
and simulation, including in the low mass region dominated by backgrounds. 

5.2.2 Results 

In total 3684 events were selected in the inclusive sample and 1595 were selected 
in the non-radiative class from data at the 183-207 GeV energy points. For the non- 




























24 


radiative sample the average efficiency in the hducial volume, taking into account all 
corrections, was 93.25 ± 0.48%, and the background from non e’''e“ —collisions 
was estimated to increase from (10 ± 2)fb at ^/s ~ 183 GeV to (13 ± 3)fb at ^/s ~ 
207 GeV, which is at most approximately 0.6% of the signal cross-section. In the inclusive 
sample the total background amounted to, at most, 9% of the signal cross-section. The 
dominant backgrounds came from four-fermion events which were misidentihed as two- 
fermion events. Gosmic ray events were estimated to account for between 0.31% and 
0.45% of the selected events. The efficiency, background etc. at each energy is given 
in Table 9. The selected events were used to measure the total cross-sections, forward- 
backward asymmetries and the differential cross-sections for non-radiative events. 


Cross-sections 


The total cross-section for the non-radiative and inclusive samples of e’''e 
events was computed from 


jj, 


a 




h47T 


eC 


( 6 ) 


where Ngei is the numbers of events selected in each sample, Nbg is the predicted number of 
background events selected in the sample for the luminosity C. The efficiency for selecting 
e'''e“ —> events in the sample is e, which includes corrections for the efficiency 

to reconstruct and identify muons and the trigger efficiency, obtained from the data. 
/ (applicable to the non-radiative cross-section only) is a factor for the migration of 
e'''e“ —events into the sample from lower centre-of-mass energy. 774 ,^ is a correction 
factor required to extrapolate the measurements from the polar angular acceptance to 
the full 47r acceptance, determined from simulated events. The measured cross-sections 
are given in Table 10 and shown in Figure 10. 

For the 207 GeV energy point separate measurements were made for the periods 
before and after the failure of the sector of the TPG. These measured cross-sections 
were combined using the BLUE technique [29] which performs a weighted average of the 
measurements taking into account correlated errors. 


Forward-backward asymmetries 


The forward-backward asymmetry for each sample was computed from 


_ Nf - CbNb 
Nf -F CbNb 


( 7 ) 


where Nf and Nb are the numbers of e’''e“ —events in which the is in the 
forward and backward hemispheres dehned as cos 6 > 0 and cos 6 < 0 respectively. 
Differences in the efficiencies for selecting events and the migration of events into the 
non-radiative samples in the forward and backward hemispheres are accounted for by 
the factor Cb, which was determined for each sample from simulated events. These dif¬ 
ferences arise from the different fractions of non-radiative and radiative return events 
in the forward and backward hemispheres reflecting the very different forward-backward 
asymmetries of these subsamples. Nf and Nb are computed from 

Nm = >)2, (VL- - . (8) 

where N^J+^- and Nlg’^ are the numbers of events selected and the predicted backgrounds 
in each hemisphere. 772,7 is a factor to extrapolate the observed number of events in 
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Energy point (GeV) 

183 

189 

192 

196 

200 

202 

205 

207 

Energy (GeV) 

182.65 

188.56 

191.60 

195.53 

199.53 

201.65 

204.85 

206.55 

Lumi (pb“^) 

52.54 

156.38 

25.79 

73.98 

83.14 

40.51 

75.55 

137.07 

No. Events 

379 

991 

167 

389 

506 

205 

373 

674 

Efficiency (%) 

92.7 

93.7 

93.1 

93.5 

93.2 

93.2 

93.6 

92.7 

Background (%) 

0.7 

0.9 

0.8 

0.8 

0.8 

0.9 

1.1 

1.0 


Table 9: Details of LEP II analysis for the channel. The table shows the 

actual centre-of-mass energy and luminosity analysed at each energy point, the number 
of events selected in the inclusive analysis and the efficiencies for selecting events in the 
non-radiative samples and the backgrounds selected in the non-radiative samples. 


e''"e ^ Mu 

(GeV) 

> 75 (GeV) 

GT/Vs > 0.85 

(Pb) 

^FB 

(pb) 

am 

^FB 

130 

24.30 ± 3.20 ± 0.83 
(20.04) 

0.450 ±0.120 ±0.002 

(0.339) 

9.70 ± 1.90 ±0.36 
(8.11) 

0.670 ±0.150 ±0.003 

(0.719) 

136 

17.00 ±2.60 ±0.58 

(17.07) 

0.560 ±0.130 ±0.002 

(0.340) 

6.60 ± 1.60 ±0.25 

(7.00) 

0.740 ±0.160 ±0.003 

(0.699) 

161 

9.30 ±1.10 ±0.32 
(10.32) 

0.390 ±0.110 ±0.002 

(0.335) 

3.60 ±0.70 ±0.13 
(4.43) 

0.430 ±0.160 ±0.003 

(0.629) 

172 

8.90 ±1.10 ±0.31 

(8.74) 

0.550 ±0.100 ±0.002 

(0.332) 

3.60 ±0.70 ±0.12 

(3.79) 

0.940 ±0.140 ±0.003 

(0.610) 

183 

8.28 ±0.46 ±0.07 

(7.56) 

0.299 ± 0.051 ± 0.003 

(0.329) 

3.61 ±0.28 ±0.03 

(3.30) 

0.588 ±0.064 ±0.001 

(0.596) 

189 

7.17 ±0.24 ±0.06 

(7.01) 

0.336 ± 0.032 ± 0.003 

(0.327) 

3.07 ±0.15 ±0.02 

(3.07) 

0.600 ±0.039 ±0.001 

(0.589) 

192 

7.37 ±0.61 ±0.07 

(6.76) 

0.275 ± 0.078 ± 0.003 

(0.327) 

2.82 ± 0.36 ± 0.02 

(2.97) 

0.636 ±0.098 ±0.001 

(0.586) 

196 

5.89 ±0.32 ±0.06 

(6.45) 

0.297 ±0.052 ±0.004 

(0.326) 

2.76 ±0.21 ±0.02 

(2.84) 

0.586 ±0.061 ±0.001 

(0.582) 

200 

6.95 ±0.33 ±0.07 

(6.16) 

0.334 ± 0.044 ± 0.003 
(0.324) 

3.08 ±0.21 ±0.02 
(2.71) 

0.548 ±0.056 ±0.001 

(0.578) 

202 

5.70 ±0.44 ±0.06 

(6.02) 

0.365 ± 0.070 ± 0.004 
(0.324) 

2.46 ±0.27 ±0.02 

(2.65) 

0.544 ±0.090 ±0.001 

(0.577) 

205 

5.46 ±0.31 ±0.06 

(5.81) 

0.304 ± 0.053 ± 0.004 

(0.323) 

2.35 ±0.19 ±0.02 

(2.56) 

0.642 ±0.061 ±0.001 

(0.574) 

207 

5.49 ±0.23 ±0.06 

(5.70) 

0.373 ± 0.039 ± 0.006 

(0.323) 

2.47 ±0.15 ±0.02 

(2.52) 

0.558 ±0.048 ±0.001 

(0.573) 


Table 10: Measured cross-sections and forward-backward asymmetries for inclusive and 
non-radiative e+e“ — events. The first uncertainty is statistical, the second sys¬ 
tematic. Numbers in brackets are the theoretical predictions of ZFITTER, which are 
estimated to have a precision o/±0.4% on and ±0.004 on Ap^. 
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the angular acceptance to 0 < |cos0| < 1. account for the expected numbers of 

events in each hemisphere with misidentihed charges, such that an event in which the 

was produced in the forward hemisphere is reconstructed as an event with the ijT in 
the backward hemisphere and vice versa. There were approximately 0.5% of events in 
which the charge of one of the particles was misidentihed, so that both particles were 
reconstructed with the same charge. In these cases, the particle with the larger relative 
uncertainty on its momentum was considered to have its charge mismeasured and its 
charge was inverted before the event was assigned to either the forward or backward 
hemisphere. From simulation this was found to identify the wrongly measured charge in 
93±5% of cases. The rate of events in which both particles had their charges mismeasured 
was found to be negligible. The measured forward-backward asymmetries are given in 
Table 10 and shown in Figure 11. 

For the 207 GeV energy point, separate measurements of the numbers of events in 
the forward and backward hemispheres were made for the periods before and after the 
failure of the sector of the TPC. These measurements were added together and an average 
forward-backward asymmetry was computed. 


Differential cross-sections 


The differential cross-sections, da/d cos 6* for the non-radiative samples of 
e’''e“ — jj,'^jj,~ events were computed in 10 bins, i, of cos 6* in the range —0.970 < cos 6* < 
0.970 using 


da/d cos 6*^ 



- % + iVi) 




(9) 


where are the number of observed events and the expected number of 

background events in each bin of cos 9, where 9 is the polar angle of the negatively charged 
muon with respect to the beam axis, is the efficiency in each bin of cos9, which 
was corrected, in all bins, for the global track reconstruction and muon identihcation 
correction factor determined for the total cross-sections and the trigger efficiency. The 
bin width in cos6' is A*. The migration factor between non-radiative and radiative return 
events, /, was computed individually for each bin from simulated events 

Migrations between neighbouring bins of cos 9 were found to be negligible. The ex¬ 
pected number of events with mismeasured charges in each bin was accounted for by 
Nl. 

Measurements in the outer bins of cos 9 were not corrected to the full angular accep¬ 
tance. Corrections were made for the experimental cuts on the polar angle of both muons 
by the factor rj which was computed from simulated events and found to be signihcant 
only in the bins with largest | cos 6^1. The measured differential cross-sections are given in 
Table 11. 

For the 207 GeV energy point separate measurements of the differential cross-sections 
in each bin were made for the periods before and after the failure of the sector of the 
TPC. The measurements were combined using the BLUE technique, taking into account 
correlated systematic errors. 


Systematic errors 

For measurements of the non-radiative cross-sections the largest single systematic 
error comes from the determination of the selection efficiency which is dominated by 
the correction factor for muon identihcation and track reconstruction efficiencies. The 
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./s ~ 189 


~ 192 

COS 6 

dcr/dcos 6 (pb) 

COS 6 

dcr/dcos^ (pb) 


Measurement 


Measurement 

[-0.97,-0.80] 

0.46 

0.41±0.13(0.14)±0.01 

1 

[-0.97,-0.80] 

0.45 

-0.01±0.00{0.33)±0.00 

[-0.80,-0.60] 

0.46 

0.47±0.13(0.12)±0.01 


[-0.80,-0.60] 

0.45 

0.21±0.21(0.31)±0.00 

[-0.60,-0.40] 

0.54 

0.54±0.14(0.14)±0.01 


[-0.60,-0.40] 

0.53 

0.83±0.42(0.33)±0.01 

[-0.40,-0.20] 

0.71 

0.43±0.12(0.16)±0.01 


[-0.40,-0.20] 

0.69 

0.42±0.30{0.38)±0.01 

[-0.20, 0.00] 

0.98 

1.21±0.21(0.19)±0.02 


[-0.20, 0.00] 

0.94 

1.09±0.49(0.45)±0.01 

[ 0.00, 0.20] 

1.33 

1.25±0.21(0.22)±0.02 


[ 0.00, 0.20] 

1.28 

0.67±0.40{0.54)±0.01 

[ 0.20, 0.40] 

1.78 

2.03±0.26(0.24)±0.02 


[ 0.20, 0.40] 

1.71 

2.98±0.80{0.60)±0.02 

[ 0.40, 0.60] 

2.33 

2.19±0.28(0.28)±0.02 


[ 0.40, 0.60] 

2.24 

2.27±0.69(0.68)±0.02 

[ 0.60, 0.80] 

2.98 

2.86±0.32(0.32)±0.03 


[ 0.60, 0.80] 

2.87 

2.30±0.70{0.77)±0.02 

[ 0.80, 0.97] 

3.70 

3.90±0.40(0.39)±0.04 


[ 0.80, 0.97] 

3.57 

3.22±0.90{0.94)±0.03 


V* ~ 183 


d(j/dcos0 (pb) 

COS 9 


Measurement 

[-0.97,-0.80] 

0.47 

-0.01±0.00(0.24)±0.00 

O 

o 

o' 

00 

o 

0.48 

0.41±0.21(0.22)±0.01 

[-0.60,-0.40] 

0.57 

1.00±0.32(0.24)±0.01 

[-0.40,-0.20] 

0.76 

1.06±0.33(0.27)±0.01 

[-0.20, 0.00] 

1.05 

1.14±0.35(0.33)±0.01 

[ 0.00, 0.20] 

1.44 

1.69±0.43(0.39)±0.02 

[ 0.20, 0.40] 

1.92 

1.79±0.43(0.44)±0.02 

[ 0.40, 0.60] 

2.52 

3.34±0.59(0.51)±0.03 

[ 0.60, 0.80] 

3.22 

3.30±0.60(0.58)±0.03 

[ 0.80, 0.97] 

4.00 

4.07±0.72(0.71)±0.04 


y/s 200 


y/s 202 

cos 6 

d(j/d COS 6 (pb) 

cos 6 

dcr/dcos^ (pb) 

SM 

Measurement 

SM 

Measurement 

[-0.97,-0.80] 

0.43 

0.53±0.20(0.18)±0.01 


[-0.97,-0.80] 

0.43 

1.43±0.48(0.26)±0.02 

[-0.80,-0.60] 

0.43 

0.26±0.13(0.17)±0.00 


[-0.80,-0.60] 

0.42 

0.54±0.27(0.24)±0.01 

[-0.60,-0.40] 

0.49 

0.83±0.23(0.18)±0.01 


[-0.60,-0.40] 

0.49 

0.25±0.18(0.25)±0.00 

[-0.40,-0.20] 

0.64 

0.71±0.22(0.20)±0.01 


[-0.40,-0.20] 

0.62 

0.26±0.19(0.29)±0.01 

[-0.20, 0.00] 

0.86 

1.10±0.28(0.24)±0.01 


[-0.20, 0.00] 

0.84 

0.56±0.28(0.34)±0.01 

[ 0.00, 0.20] 

1.17 

1.44±0.32(0.28)±0.01 


[ 0.00, 0.20] 

1.14 

0.71±0.32(0.40)±0.01 

[ 0.20, 0.40] 

1.56 

1.66±0.33(0.32)±0.02 


[ 0.20, 0.40] 

1.52 

2.45±0.58(0.45)±0.02 

[ 0.40, 0.60] 

2.04 

2.72±0.42(0.36)±0.02 


[ 0.40, 0.60] 

1.99 

1.45±0.44(0.51)±0.02 

[ 0.60, 0.80] 

2.61 

2.80±0.43(0.41)±0.02 


[ 0.60, 0.80] 

2.55 

1.87±0.50(0.58)±0.02 

[ 0.80, 0.97] 

3.25 

3.18±0.50(0.50)±0.03 


[ 0.80, 0.97] 

3.18 

3.14±0.71(0.71)±0.03 


y/s ^ 196 


d(j/d COS 9 (pb) 

cos 9 

SM 

Measurement 

[-0.97,-0.80] 

0.44 

0.68±0.24(0.19)±0.01 

O 

o 

o' 

00 

o 

0.44 

0.51±0.19(0.18)±0.01 

o 

o 

o'' 

o 

o 

0.51 

0.65±0.22(0.19)±0.01 

[-0.40,-0.20] 

0.66 

0.22±0.13(0.22)±0.01 

[-0.20, 0.00] 

0.90 

0.84±0.26(0.26)±0.01 

[ 0.00, 0.20] 

1.22 

1.71±0.37(0.31)±0.02 

[ 0.20, 0.40] 

1.63 

1.33±0.32(0.35)±0.01 

[ 0.40, 0.60] 

2.14 

2.08±0.39(0.39)±0.02 

[ 0.60, 0.80] 

2.74 

2.92±0.46(0.45)±0.02 

[ 0.80, 0.97] 

3.41 

2.75±0.49(0.54)±0.03 


y* ~ 205 


^ ~ 207 

COS 9 

d(j/d COS 9 (pb) 

COS 9 

dcr/dcos^ (pb) 

SM 

Measurement 

SM 

Measurement 

[-0.97,-0.80] 

0.42 

0.67±0.24(0.19)±0.01 


[-0.97,-0.80] 

0.42 

0.47±0.15(0.14)±0.01 

[-0.80,-0.60] 

0.42 

0.35±0.16(0.17)±0.00 


[-0.80,-0.60] 

0.41 

0.47±0.13(0.13)±0.00 

[-0.60,-0.40] 

0.47 

0.20±0.12(0.18)±0.00 


[-0.60,-0.40] 

0.47 

0.42±0.13(0.13)±0.00 

[-0.40,-0.20] 

0.60 

0.34±0.16(0.20)±0.01 


[-0.40,-0.20] 

0.59 

0.54±0.14(0.15)±0.01 

[-0.20, 0.00] 

0.81 

0.57±0.20(0.24)±0.01 


[-0.20, 0.00] 

0.80 

0.71±0.16(0.18)±0.01 

[ 0.00, 0.20] 

1.10 

0.95±0.27(0.28)±0.01 


[ 0.00, 0.20] 

1.08 

1.00±0.20(0.21)±0.01 

[ 0.20, 0.40] 

1.47 

1.33±0.31(0.32)±0.01 


[ 0.20, 0.40] 

1.44 

1.39±0.23(0.24)±0.01 

[ 0.40, 0.60] 

1.92 

1.40±0.32(0.37)±0.02 


[ 0.40, 0.60] 

1.88 

2.61±0.32(0.27)±0.02 

[ 0.60, 0.80] 

2.46 

2.16±0.40(0.42)±0.02 


[ 0.60, 0.80] 

2.41 

1.93±0.28(0.31)±0.02 

[ 0.80, 0.97] 

3.06 

4.23±0.61(0.51)±0.03 


[ 0.80, 0.97] 

3.01 

2.65±0.36(0.38)±0.02 


Table 11: Differential cross-sections for non-radiative events at cen- 

tre-of-mass energies from 183 to 207 GeV. The tables show the bins, the predictions 
of the Standard Model and the measurements. The errors quoted are the statistical and 
experimental systematic errors. The statistical errors are shown as the measured errors 
and, in brackets, the expected errors, computed from the square root of the observed and 
expected numbers of events respectively. 
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e+e —p 




Source 

Aa/a 

(non-rad.) 

Aa/a 

(inclus.) 

AAfb 

(non-rad.) 

AAfb 

(inclus.) 

Efficiency 

51 

51 

9 

20 

\/p reconstruction 

14 

20 

2 

0 

Charge misidentification 

- 

- 

6 

4 

QED 

2 

9 

1 

0 

Backgrounds 

12 

70 

6 

51 

Luminosity 

55 

55 

- 

- 

Total uncertainty 

77 

105 

12 

55 

Correlated 

75 

103 

11 

53 

Uncorrelated 

16 

17 

5 

17 


Table 12: Systematic uncertainties on the measured cross-sections and forward-backward 
asymmetries for data taken at ~ 207 GeV. All numbers in units of 10~^. Separate 
measurements were made before and after the TPC sector failure, and then combined. 
There is thus no explicit systematic uncertainty related to the TPC sector instability. 
The correlated error component includes errors correlated between energies and channels 
and those correlated with other LEP experiments. 
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relative error on the correction factor is 0.49% and is correlated between all energies and 
between all bins in the estimation of the differential cross-sections. Additional, smaller 
sources of error on the efficiency come from the statistics of simulated event samples and 
the estimated trigger efficiency. 

The uncertainty on the ISR migration factors arising from misreconstruction of 
has two components: an uncorrelated component from the statistics of simulated 
e'''e“ —events; and an uncertainty estimated from switching off the simulation 
of the beam energy spread in the simulated sample, which was found by comparing the 
change in the migration factor from the standard analysis. This was determined for a 
single energy and applied as a correlated error over all energies. These uncertainties were 
estimated bin-by-bin for the differential cross-sections. 

The uncertainty on the angular acceptance corrections comes from the statistics of the 
simulated e’''e“ — /jT samples. The uncertainty on the forward-backward asymmetry 
and differential cross-sections due to charge misidentihcation comes from the estimation 
of the efficiency for correctly identifying which particle had the mismeasured charge in 
such events. The uncertainty amounts to, at most, 0.0006 on the non-radiative forward- 
backward asymmetries. This introduces an anticorrelation between bins in the forward 
and backward hemispheres in the differential cross-sections. 

The uncertainties arising from missing higher orders in the QED radiative corrections 
in the simulation were estimated by taking half the difference in the measurements ob¬ 
tained using KK simulated events in which the highest order corrections available in the 
program were either included or excluded. 

Uncertainties on the backgrounds come from the statistics of the simulated samples and 
from the number of cosmic rays failing vertex-selection cuts used to estimate the residual 
cosmic ray background, which are uncorrelated between energy points. The hnal source 
comes from the theoretical precision on the total cross-section of the different simulated 
background samples, which is correlated between energy points and is signihcant for the 
measurement of the inclusive cross-sections and forward-backward asymmetries. 

Uncertainties from the evaluation of the luminosity were included in the measurements 
of the cross-sections. The systematic errors are summarised in Table 12. 

5.3 T+T” final states 

The analysis of tan pair production in 1997-98 LEP runs was presented in [4]. Here 
we present a complete reanalysis of those data together with the new analysis of data 
taken in 1999-2000. Several aspects of the analysis have changed with respect to the 
previous publication; improved track reconstruction algorithm, optimised event selection 
procedure and a better description of signal and background processes with new event 
generators. All the data taken in the 1997-2000 runs were analysed with a homogeneous 
procedure which is described below. The data taken above the Z resonance in the 1995-96 
runs [3] were not reanalysed since the improvements of the new analysis are negligible 
compared to the large statistical uncertainty of those measurements. 

5.3.1 Analysis 
Run selection 

Runs in which detector components essential for the analysis were not fully operational 
were rejected to ensure a high quality of the data used in the analysis. This included the 
requirement that the TPC, HPC and Forward Electromagnetic Calorimeters (FEMC) 
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had a performance close to optimal. In addition, from the two pairs of detectors VD and 
Inner Detector (ID), and the forward tracking chambers, FCA and FCB, at least one of 
the two detectors in both pairs had to be operational. About 2.5% of total integrated 
luminosity was rejected by this procedure. 

Track and photon selection 

The r-pair selection was based on the reconstructed kinematic properties of the events. 
Therefore it was important that only well reconstructed tracks and photons were used in 
the analysis. To ensure this the following selection was applied to the tracks found by 
the reconstruction program: 

• The particle momentum had to be at least 100 MeV/c; 

• The error of the momentum reconstruction had to be less than the absolute value 
of the momentum; 

• The track extrapolation had to originate from the region around the nominal inter¬ 
action point (less than 5 cm in the plane perpendicular to the beam, less than 10 
cm along the z axis); 

• The track had to be seen either in VD or in ID; 

• The track had to be seen in at least one tracking device beyond VD and ID 
(namely TPC, OD, FCA, FCB). 

The last condition was not applied to the tracks close to the azimuthal boundaries of 
the TPC sectors, nor was it applied to the tracks pointing to the broken TPC sector for 
the last period of 2000 data taking. In these cases, if the last condition was not fulhlled 
the track was still accepted if it had a measured RICH tracking point or if there was any 
deposition in an electromagnetic or hadron calorimeter within 10 degrees from the track 
direction. 

An unassociated electromagnetic cluster was accepted as a photon if it had more than 
0.5 GeV of energy deposition and was found more than one degree from the nearest track. 
Identihed electron-positron pairs compatible with a photon conversion were considered 
as photons rather than charged particles. 

The event selection was based on the tracks and photons accepted by these procedures. 

Event selection 

The selection of r-pairs was similar to the one presented in [4]. Several cuts were 
changed to optimise the product of efficiency and purity. 

Events with fewer than 2 tracks were rejected. The charged particle tracks were 
grouped into two jets using the LUCLUS algorithm [20]. Tracks in each jet were con¬ 
sidered as decay products of a r-candidate. The most energetic charged particle in each 
jet was called the leading charged particle and its track was called the leading track. The 
momentum direction of the r-candidate was approximated by the vector sum of momenta 
of its charged decay products. 

Under the assumption that a single ISR photon is emitted along the beam direction, 
the r-candidate momenta can be estimated (neglecting the r mass) from: 

o r ^ 2__ P 

* sin 9i + sin 62 + \ sin ( 6^1 -|- 6 ^ 2 ) | 


( 10 ) 
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where 6 i are the polar angles of the tau-candidate momenta, and j denotes the r candidate 
opposite to the candidate i. The reduced centre-of-mass energy, V^, is then given by: 

s' _ sin 9i + sin 62 - \ sin {9i + 6 ^ 2 ) | , , 

s sin 9i + sin 6^2 + | sin ( 6*1 + 6 * 2 ) | 

The r pairs were selected using the set of cuts listed below: 

• To reject multi-hadronic events the number of selected tracks Nch had to be less 
than seven: N^h < 7. For events in which both jets contained more than one track 
the invariant mass of each jet (computed using the charged particles, which were 
assumed to be pious) had to be less than 2.5 GeV/c^; 

• The background from WW pairs, ZZ pairs and the e’''e“ —> e’''e“ff process was 
suppressed by requiring that the acoplanarity of the two leading tracks had to be 
less than 12 degrees: 9acop < 12 °; 

• To suppress background from two-photon collisions the visible energy (computed as 

a scalar sum of charged particle momenta and photon energies) had to be larger 
than 20% of the nominal centre-of-mass energy: Eyis > 0.2^. The cut value was 
increased to 0.225y^ for e’''e“ and candidates (events with exactly two tracks 

both satisfying loose electron identihcation or muon identihcation criteria); 

• Additional suppression of two-photon background was achieved using a cut on the 

transverse momentum (computed as the projection onto the r -0 plane of the 
vector sum of all charged particle momenta): Pr < 2.5 GeV/c. The cut value was 
increased to 8 GeV/c for e’''e“ and candidates; 

• The background from e’''e“ —/i■*'/!“ and partially from e^e“ —e’''e“ events was 

rejected using the “radial momentum” variable, Prad = ^(pi/Pi)^ + (P 2 /P 2 )^) 
where Pi are the momenta of leading charged particles and P[ are the r momenta 
estimated using eqn. 10 . The selection cut was Prad < 0.95 for candidates 

and Prad < 1-2 for all other events; 

• The bulk of the e+e“ —> e+e“ background was rejected using the depositions in 
the electromagnetic calorimeters HPG and FEMG. The “reduced electromagnetic 
energies” e* were dehned for each jet as e* = Ei/EJ, where Ei are the total elec¬ 
tromagnetic energy depositions in cones of half angle 30° around the two leading 
charged particles and EJ are the estimated r energies, computed from the estimated 
momenta P[ of eqn. 10. For the events contained within the HPG acceptance both 
reduced electromagnetic energies had to satisfy the condition < 0.85. In addition, 
the distance between the point with coordinates (ci, 62 ) and the point ( 1 , 1 ) expected 
for e+e“ —e^e“ events had to be more than 0.5. In the forward part of the detec¬ 
tor the e+e“ —e+e“ background was significantly higher than in the barrel part. 
Therefore the selection cuts were slightly tighter for the events contained within the 
FEMG acceptance and an additional cut was applied: the energy deposited beyond 
the hrst layer of the HGAL had to be non-zero; 

• The background from cosmic-ray events was rejected using the impact parameters 

of the two leading tracks in the T-(j) plane and along the z axis The 

values of had to be not simultaneously larger than -|-3 cm or simultaneously 

smaller than -3 cm and at least one of or R^^^ had to be less than 3 mm in 

absolute value. In addition, for the two-track events, the difference \Zl^^ — 

had to be less than 3 cm; 

• The remaining background from e^e“ —> e^e“, e^e“ —and cosmic-ray events 
was reduced by rejecting back-to-back events with very low acollinearity between the 
two leading tracks: 9acol < 0.3°. 
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Figure 5: The reconstructed a/s'/a/s distribution (left) and the reconstructed a/s' distri¬ 
bution (right) for e+e“ —> r+r“ events. Both plots show data from 1997-2000. 

After the selection procedure the residual background level was relatively low (10- 
20%) for non-radiative events and for radiative return events. However the region of 
intermediate values of reconstructed a/s' was dominated by background from 77 —>• e’''e“ 
events. Therefore all e+e“ candidate events were rejected in the region \Jl/3 < /^/s < 

0.85. 

Events were associated to the non-radiative or inclusive samples according to the value 
of estimated using eqn. 11. The distribution of the reduced centre-of-mass energy, 
for selected r pair events from centre-of-mass energies of 183 to 207 GeV, is shown in 
Figure 5. The ratio v^/a/s is shown for the non-radiative sample, while for the events 
with a hard ISR photon the value of a/s', which approximates the invariant mass of the 
hnal state leptons, is shown. 

Estimation of the selection efficiency and residual background 

Unlike in the cases of e’''e“ —^ e’''e“ and e’''e“ —>■ pi~, the selection efficiency of 
e'''e“ —> T^T~ events cannot be estimated or verihed from real data. The estimation 
of the selection efficiency is entirely based on simulations. To ensure good agreement 
between real and simulated data an extensive study was performed using test samples of 
real and simulated events. Where necessary, the simulation was corrected by introducing 
calibration constants, smearing distributions, etc. 

The (^-dependent) energy scale and energy resolution of the electromagnetic calorime¬ 
ters were calibrated using well reconstructed e’'“e“ —>■ e’''e“ events selected from high en¬ 
ergy runs. Non-linearities were checked using e'''e“ —> e'''e“ events selected from Z runs 
and 77 —> e’''e“ events at high energies. The corresponding re-scaling and smearing were 
applied to electromagnetic energies in simulated events. Also, a small forward-backward 
asymmetry of the electromagnetic calorimeters was found and corrected for. 
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The absolute scale and resolution of charged-particle momentum measurements were 
calibrated using e’''e“ —/i■*■/!“ events from Z runs. The real data needed a small 6 - 
dependent correction to the measured 1 /pr, while simulated data needed a smearing 
of the 1/pt distribution. The effect of this calibration was checked with non-radiative 
e“''e“ —events from high energy data. A rather small (about 2a) discrepancy was 
found and taken as a systematic uncertainty. 

The momentum dependence of the muon chamber efficiency was studied using test 
samples of muons selected with very tight HCAL criteria. The e’''e“ —pT events 
from Z and high energy runs were used to select muons with ~45 GeV/c and ~100 
GeV/c momenta, while 77 —and r —/rz/z/ events from high energy runs were 
used for momentum regions of 0-45 and 45-100 GeV/c. Both the mean efficiency and 
momentum dependence were corrected in the simulation. 

Other calibrations included the electron punch-through in the hrst layer of the HGAL, 
the specihc energy loss (dE/dx) measurements from the TPG and the LEP energy spread. 
After applying all necessary corrections and calibrations the selection efficiency was cal¬ 
culated from simulated data as the fraction of events generated above the nominal 
cut which pass the selection criteria. 

The residual background level was normalised using the real data. For each type of 
background a selection cut most sensitive to this background source was chosen and all 
other cuts were applied to the real and simulated data. The events failing the “sensitive 
cut” were then used to normalise the background level predicted by simulation to the 
real data. As an example the distributions of acoplanarity (sensitive to four-fermion 
background) and radial momentum (sensitive to e’''e“ — p~) are presented in Figure 6 . 
The full statistics of the 1997-2000 runs are shown. Arrows show the selection cuts; events 
to the left of the arrows were selected as r candidates and events to the right were used 
for background level normalisation. 

Due to this normalisation procedure the simulation is used only to predict the shapes 
of distributions, while the absolute background level is estimated from the real data. The 
size of the corrections applied to the simulation was of the order of 5-15%. The estimated 
residual backgrounds from different channels are presented in Table 13 for non-radiative 
and inclusive samples selected from 2000 data. For other years the background levels 
were similar. The estimation of background uncertainties is discussed in section 5.3.2. 

The feed-through of e+e“ —r+r“ events from lower values of V ^/\/s into the non- 
radiative sample was estimated from simulated events to be approximately 4.4%, almost 
independent of ^/s. 

5.3.2 Results 


Cross-sections 

The total cross-sections of r pair production for the non-radiative and inclusive samples 
were calculated as follows: 


(Ttt 


iNsel-mg)il-f) 

eC 


( 12 ) 


Here C is the total integrated luminosity for the corresponding energy point, Nsei is 
the number of events selected in the corresponding sample, / (applicable to the non- 
radiative cross-section only) is the fraction of feed-through events estimated from the r- 
pair simulation and e is the selection efficiency, computed within the full 47 r acceptance. 
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Figure 6: The acoplanarity distribution (left) and the distribution of the variable Prad 
(right). The acoplanarity distribution is sensitive to the four-fermion backgrounds while 
the Prad distribution is sensitive to the yi~ background. Both plots show data 

from 1997-2000. The dominant backgrounds have been normalised to obtain good agree¬ 
ment between simulation and data in the region dominated by background. 


e e 


T^T~ 


Background 


Non-radiative 


Inclusive 


source 

e'*'e~ —*• e^e~ 
e’''e“ — 
e’''e“ —> qq 

77 —>■ T~^T~ 

77 —>• e’''e“ 

77 ^ /i+/i“ 

77 ^ qq 

e+e- ^ W+W- 
e+e- ^ ZZ 
e+e“ —»• e+e“ff 
Cosmic rays 
Total 


background 


4.50 

± 

0.54 

1.05 

± 

0.15 

0.19 

± 

0.03 

0.02 

± 

0.01 

0.11 

± 

0.05 

0.07 

± 

0.04 

0.00 

± 

0.04 

2.93 

± 

0.50 

1.09 

± 

0.19 

2.19 

± 

0.38 

0.28 

± 

0.22 


12.43 ± 0.89 


background 


3.37 

± 

0.37 

1.18 

± 

0.16 

1.40 

± 

0.14 

0.23 

± 

0.05 

1.89 

± 

0.31 

0.23 

± 

0.06 

0.10 

± 

0.06 

3.15 

± 

0.54 

2.04 

± 

0.35 

5.50 

± 

0.94 

0.24 

± 

0.15 


19.34 ± 1.26 


Table 13: Residual background level for 2000 data (in % relative to the number of selected 
T-pair candidates). 
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e“'"e —>■ 







Energy point (GeV) 

183 

189 

192 

196 

200 

202 

205 

207 

Energy (GeV) 

182.66 

188.63 

191.58 

195.51 

199.51 

201.64 

204.89 

206.56 

Lumi (pb“^) 

53.11 

152.67 

25.14 

75.99 

82.65 

40.40 

81.49 

136.39 

No. Events 

228 

551 

79 

232 

269 

119 

273 

401 

Efficiency (%) 

52.1 

52.2 

52.0 

52.0 

52.0 

52.0 

51.7 

51.7 

Background (%) 

12.2 

11.8 

14.6 

12.6 

13.6 

13.8 

12.3 

12.5 


Table 14: Details of LEP II analysis for e+e“ —r+r“ channel. The table shows the 
actual centre-of-mass energy and luminosity analysed at each energy point, the number 
of events selected in the inclusive analysis and the efficiencies for selecting events in the 
non-radiative samples and the backgrounds selected in the non-radiative samples. 


e"*”© — t'*”t 

±7 

(GeV) 

> 75 (GeV) 

GT/Vs > 0.85 

o-TT (pb) 

^FB 

CTtt (pb) 

^FB 

130 

22.20 ± 4.60 ± 1.56 

(20.30) 

0.310 ±0.170 ±0.020 
(0.337) 

10.20 ±3.10 ±0.72 

(8.31) 

0.730 ±0.170 ±0.020 
(0.719) 

136 

17.70 ± 3.90 ± 1.24 

(17.29) 

0.260 ±0.190 ±0.020 

(0.338) 

8.80 ± 3.00 ± 0.62 

(7.17) 

0.490 ± 0.230 ± 0.020 

(0.699) 

161 

11.70 ±1.80 ±0.82 

(10.44) 

0.390 ±0.120 ±0.020 

(0.332) 

5.10 ± 1.20 ±0.36 
(4.54) 

0.920 ± 0.080 ± 0.020 

(0.628) 

172 

11.20 ±1.80 ±0.79 

(8.83) 

0.190 ±0.140 ±0.020 
(0.329) 

4.50 ± 1.10 ±0.32 

(3.89) 

0.130 ±0.200 ±0.020 
(0.610) 

183 

8.73 ± 0.70 ± 0.20 

(7.64) 

0.400 ± 0.074 ± 0.012 

(0.326) 

3.29 ±0.38 ±0.07 

(3.39) 

0.671 ±0.080 ±0.012 

(0.596) 

189 

7.23 ±0.38 ±0.17 

(7.08) 

0.448 ±0.047 ±0.012 

(0.324) 

3.11 ±0.22 ±0.07 

(3.15) 

0.697 ±0.048 ±0.011 

(0.589) 

192 

6.16 ±0.89 ±0.15 

(6.83) 

0.415 ±0.134 ±0.012 

(0.324) 

2.50 ±0.48 ±0.06 

(3.04) 

0.578 ±0.150 ±0.011 

(0.586) 

196 

5.97 ±0.51 ±0.14 
(6.52) 

0.156 ±0.080 ±0.012 
(0.323) 

2.89 ± 0.30 ± 0.06 
(2.91) 

0.465 ±0.083 ±0.011 
(0.582) 

200 

6.50 ±0.50 ±0.15 
(6.22) 

0.226 ± 0.073 ± 0.012 
(0.322) 

2.61 ±0.27 ±0.06 

(2.78) 

0.540 ± 0.080 ± 0.011 

(0.578) 

202 

5.74 ±0.68 ±0.14 

(6.08) 

0.314 ±0.110 ±0.012 

(0.321) 

2.55 ±0.38 ±0.06 
(2.72) 

0.464 ±0.122 ±0.011 

(0.576) 

205 

6.94 ±0.52 ±0.16 

(5.86) 

0.434 ± 0.070 ± 0.012 

(0.320) 

2.80 ±0.28 ±0.06 
(2.62) 

0.709 ±0.068 ±0.011 

(0.574) 

207 

5.95 ±0.38 ±0.14 
(5.76) 

0.347 ± 0.060 ± 0.012 
(0.320) 

2.53 ±0.21 ±0.06 
(2.58) 

0.666 ±0.059 ±0.011 
(0.572) 


Table 15: Measured cross-sections and forward-backward asymmetries for inclusive and 
non-radiative e’''e“ — events. The first uncertainty is statistical, the second sys¬ 
tematic. Numbers in brackets are the theoretical predictions of ZFITTER, which are 
estimated to have a precision of ±0.4% on cTt-t and ±0.004 on Apg. 
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which takes into account the effect of selection cuts (including the cut on \/^) and 
the trigger efficiency. The trigger efficiency was estimated from the real data using the 
redundancy of the DELPHI trigger system. Its typical value was about 99.9% with 
typical uncertainty of about 0.05%. For the 2000 data taken with one TPC sector broken 
the trigger efficiency value was 98.9±0.9% (for the whole detector). The values of the 
selection efficiency are summarised in Table 14. The residual number of background 
events was estimated from simulation and normalised to the real data as discussed in 
the previous section. The measured cross-sections are presented in Table 15 and shown 
in Figure 10. 


Forward-backward asymmetries 

The forward-backward asymmetry of r pair production was calculated as 


A 


TT _ 

FB — 


ap — cb 
O'f + O'b 


(13) 


The forward cross-section ap (0° < 9^-- < 90°) and the backward cross-section as 
(90° < 6r- < 180°) were determined in a similar way to the determination of the total 
cross-section arr- The r charge and the direction of the r momentum were estimated 
from the charges and momenta of its charged decay products. To reduce charge misiden- 
tification, only events with at least one jet consisting of exactly one charged particle [l-N 
topology events) were used in this analysis. The charge of that single particle was assigned 
to the corresponding r and an opposite charge was assigned to the other r regardless of 
the charge of the multi-track jet. The 1-1 topology events with two reconstructed tracks 
of the same charge {like-sign events) were excluded from the analysis. For the period of 
broken TPC sector in 2000 the tracks detected only in the VD and ID were not used for 
the charge identification. 

Misidentification of the r charge results in a reduction of the absolute value of ApB- 
The probability of charge misidentification was estimated from the real data using the rate 
of the like-sign events. The corresponding correction io Afb was found to be -1-0.0068 ± 
0.0034 (accidentally the same for non-radiative and inelusive samples). The corrected 
values oi Afb are summarised in Table 15 and shown in Figure 11. 


Differential cross-sections 

In addition to the measurements of the forward-backward asymmetries, differential 
cross-sections, d(T/dcos6* for r-pair production were determined. The region of angular 
acceptance of this analysis (—0.96 < cosO < 0.96) was divided into ten bins with respect 
to the polar angle of the momentum direction of the negatively charged tan lepton. The 
cross-section for every bin was determined similarly to the determination of the total 
cross-section and the correction for the charge misidentification effect was applied. The 
measured differential cross-sections are presented in Table 16. 


Systematic errors 

A breakdown of the systematic uncertainties for cross-sections and asymmetries of 
non-radiative and inelusive samples is presented in Table 17. The numbers are given for 
the results of 2000; for other years the uncertainties were of a similar size. 

The most important sources of systematic uncertainty were the choice of track- and 
event-selection cuts, background level and detector calibration. To estimate the effect of 
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7 -+^- 


yCs ~ 183 


~ 189 

cos 6 

dcr/dcos0 (pb) 

COS 9 

dcr/dcos 6 (pb) 


Measurement 


Measurement 

'o' 

00 

o 

o 

o 

o 

0.48 

-0.17±0.00(0.50)±0.03 

1 

[-0.96,-0.80] 

0.47 

0.19±0.20(0.27)±0.03 

'o' 

o 

o 

o'" 

00 

o 

0.49 

0.54±0.37(0.36)±0.03 


[-0.80,-0.60] 

0.47 

0.13±0.14(0.21)±0.03 

[-0.60,-0.40; 

0.59 

0.42±0.28(0.32)±0.02 


[-0.60,-0.40] 

0.56 

0.35±0.16(0.18)±0.03 

[-0.40,-0.20; 

0.78 

0.44±0.28(0.36)±0.02 


[-0.40,-0.20] 

0.73 

0.74±0.21{0.21)±0.03 

[-0.20, o.oo; 

1.08 

1.22±0.48(0.46)±0.04 


[-0.20, 0.00] 

1.00 

0.78±0.24(0.26)±0.04 

[ 0.00, o.2o; 

1.47 

1.51±0.53(0.52)±0.05 


[ 0.00, 0.20] 

1.36 

1.78±0.34(0.30)±0.04 

[ 0.20, 0.40; 

1.97 

1.86±0.55(0.55)±0.06 


[ 0.20, 0.40] 

1.83 

2.31±0.35(0.31)±0.05 

[ 0.40, 0.60; 

2.58 

2.30±0.61(0.64)±0.07 


[ 0.40, 0.60] 

2.39 

2.59±0.37(0.36)±0.06 

[ 0.60, 0.80; 

3.30 

3.76±0.94(0.89)±0.12 


[ 0.60, 0.80] 

3.06 

2.61±0.48(0.51)±0.11 

[ 0.80, 0.96; 

4.08 

4.77±1.40(1.31)±0.15 


[ 0.80, 0.96] 

3.78 

3.14±0.67(0.73)±0.13 


yS ~ 192 


dcr/dcos0 (pb) 

COS 9 


Measurement 

[-0.96,-0.80] 

0.46 

0.60±0.74(0.69)±0.02 

[-0.80,-0.60] 

0.46 

0.22±0.41(0.52)±0.03 

[-0.60,-0.40] 

0.54 

0.45±0.41(0.45)±0.01 

[-0.40,-0.20] 

0.71 

0.44±0.41(0.50)±0.02 

[-0.20, 0.00] 

0.97 

0.85±0.60(0.63)±0.03 

[ 0.00, 0.20] 

1.32 

1.48±0.76(0.72)±0.03 

[ 0.20, 0.40] 

1.76 

1.77±0.77(0.77)±0.04 

[ 0.40, 0.60] 

2.30 

1.67±0.77{0.88)±0.05 

[ 0.60, 0.80] 

2.95 

2.57±1.18(1.25)±0.07 

[ 0.80, 0.96] 

3.65 

0.25±0.74(1.75)±0.09 


y/s ^ 200 


y/s ^ 202 

COS 9 

dcr/d COS 9 (pb) 

cos 9 

dcr/d COS 0 (pb) 

SM 

Measurement 

SM 

Measurement 

[-0.96,-0.80] 

0.44 

0.32±0.32(0.36)±0.01 


[-0.96,-0.80] 

0.44 

-0.13±0.00(0.52)±0.01 

[-0.80,-0.60] 

0.44 

0.81±0.35(0.28)±0.03 


[-0.80,-0.60] 

0.43 

0.83±0.50(0.39)±0.03 

[-0.60,-0.40] 

0.51 

0.50±0.23(0.24)±0.01 


[-0.60,-0.40] 

0.50 

0.42±0.31(0.34)±0.01 

[-0.40,-0.20] 

0.65 

0.49±0.24(0.27)±0.02 


[-0.40,-0.20] 

0.64 

0.60±0.36(0.38)±0.02 

[-0.20, 0.00] 

0.88 

0.98±0.35(0.33)±0.03 


[-0.20, 0.00] 

0.86 

1.54±0.61{0.47)±0.03 

[ 0.00, 0.20] 

1.20 

0.93±0.34(0.38)±0.03 


[ 0.00, 0.20] 

1.17 

0.65±0.42(0.54)±0.03 

[ 0.20, 0.40] 

1.60 

1.96±0.44(0.41)±0.04 


[ 0.20, 0.40] 

1.56 

2.46±0.70(0.58)±0.03 

[ 0.40, 0.60] 

2.09 

1.86±0.44(0.47)±0.05 


[ 0.40, 0.60] 

2.04 

2.01±0.65(0.66)±0.05 

[ 0.60, 0.80] 

2.68 

2.20±0.61(0.66)±0.07 


[ 0.60, 0.80] 

2.62 

1.11±0.69(0.93)±0.07 

[ 0.80, 0.96] 

3.32 

2.89±0.87(0.92)±0.08 


[ 0.80, 0.96] 

3.24 

2.28±1.12(1.30)±0.08 


196 


dcr/d COS 9 (pb) 

cos 9 

SM 

Measurement 

[-0.96,-0.80] 

0.45 

0.35±0.34(0.39)±0.01 

O 

o 

o' 

00 

o 

0.45 

0.36±0.27(0.29)±0.03 

[-0.60,-0.40] 

0.52 

1.02±0.33(0.25)±0.01 

[-0.40,-0.20] 

0.68 

0.73±0.29(0.28)±0.02 

[-0.20, 0.00] 

0.92 

1.31±0.41(0.35)±0.03 

[ 0.00, 0.20] 

1.25 

1.03±0.37(0.40)±0.03 

[ 0.20, 0.40] 

1.68 

1.38±0.40(0.43)±0.04 

[ 0.40, 0.60] 

2.19 

2.24±0.50(0.50)±0.05 

[ 0.60, 0.80] 

2.81 

2.85±0.71(0.70)±0.07 

[ 0.80, 0.96] 

3.48 

2.92±0.91(0.98)±0.08 


y* ~ 205 


yi ~ 207 

COS 9 

dcr/d COS 9 (pb) 

COS 9 

dcr/d COS 0 (pb) 

SM 

Measurement 

SM 

Measurement 

o 

00 

o 

o 

0.43 

-0.12±0.00(0.36)±0.01 


[-0.96,-0.80] 

0.43 

0.30±0.23(0.27)±0.01 

o 

o 

o' 

00 

o 

0.42 

0.48±0.23(0.27)±0.02 


[-0.80,-0.60] 

0.42 

0.00±0.15(0.21)±0.02 

[-0.60,-0.40] 

0.48 

0.44±0.23(0.23)±0.01 


[-0.60,-0.40] 

0.48 

0.50±0.18(0.18)±0.01 

[-0.40,-0.20] 

0.62 

0.53±0.26(0.26)±0.01 


[-0.40,-0.20] 

0.61 

0.60±0.20(0.20)±0.01 

[-0.20, 0.00] 

0.83 

0.48±0.28(0.33)±0.02 


[-0.20, 0.00] 

0.82 

0.62±0.22(0.25)±0.02 

[ 0.00, 0.20] 

1.13 

1.29±0.39(0.37)±0.02 


[ 0.00, 0.20] 

1.11 

1.66±0.32(0.28)±0.02 

[ 0.20, 0.40] 

1.50 

1.30±0.37(0.40)±0.03 


[ 0.20, 0.40] 

1.48 

1.61±0.32{0.31)±0.03 

[ 0.40, 0.60] 

1.97 

2.25±0.49(0.46)±0.04 


[ 0.40, 0.60] 

1.93 

1.54±0.32{0.35)±0.04 

[ 0.60, 0.80] 

2.52 

2.87±0.66(0.64)±0.06 


[ 0.60, 0.80] 

2.48 

1.72±0.42(0.49)±0.06 

[ 0.80, 0.96] 

3.13 

4.99±1.08(0.90)±0.07 


[ 0.80, 0.96] 

3.07 

4.32±0.81(0.69)±0.07 


Table 16: Differential cross-sections for non-radiative e'^e~ ^ events at cen- 
tre-of-mass energies from 183 to 207 GeV. The tables show the bins, the predictions 
of the Standard Model and the measurements. The errors quoted are the statistical and 
experimental systematic errors. The statistical errors are shown as the measured errors 
and, in brackets, the expected errors, computed from the square root of the observed and 
expected numbers of events respectively. 
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e+e —r^r 

Source 

A(t/(T 

(non-rad.) 

A(t/(T 

(inclus.) 

AAfb 

(non-rad.) 

AAfb 

(inclus.) 

Track Selection 

72 

35 

34 

48 

Event Selection 

127 

104 

51 

59 

Detector Calibration 

91 

86 

58 

50 

Background Level 

102 

157 

39 

43 

Light pair subtraction 

9 

12 

4 

3 

Trigger Efficiency 

9 

9 

0 

0 

y/p Reconstruction 

25 

3 

14 

1 

Eeed-through 

32 

2 

17 

1 

Charge Misidentihcation 

0 

0 

34 

34 

Simulation Statistics 

21 

19 

22 

22 

QED 

7 

14 

1 

17 

Tau Polarisation 

33 

14 

16 

23 

Tau branchings 

6 

11 

0 

0 

Luminosity 

55 

55 

- 

- 

Total uncertainty 

216 

220 

105 

113 

Correlated 

208 

213 

97 

106 

Uncorrelated 

59 

53 

41 

39 

TPC sector instability 

127 

90 

41 

16 


Table 17: Systematic uncertainties for the 2000 data. All numbers in units of 10~‘^. The 
total uncertainties do not include the error due to TPC sector instability which applies 
only to a part of the 2000 data. The correlated error component includes errors correlated 
between energies and channels and those correlated with other LEP experiments. 
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track selection every selection cut was varied in a wide range and the full analysis chain 
was repeated. If the observed change of result exceeded the expected statistical fluctua¬ 
tion (la) then the quadratic difference between the change and the expected fluctuation 
was taken as the systematic error, otherwise no systematic error was assigned. The 
largest uncertainty came from the impact parameters, , and from the combination 

of detectors which was used in the track £t. 

The influence of event-selection cuts was checked in a similar way. Every event- 
selection cut was varied typically by the experimental resolution on the corresponding 
variable. The largest contribution came from the cuts on Eyis, Prad and electromag¬ 
netic energy. To reduce the effect of statistical fluctuations the full statistics of 1997-2000 
were used to estimate the uncertainty from track and event selections and therefore these 
systematic errors were common and correlated for the data of different years. 

The uncertainty of the detector calibration was mainly due to the limited statistics of 
test samples selected from real data of Z and high energy runs to estimate various correc¬ 
tions to the detector simulation. This part of the uncertainty was uncorrelated between 
different years. The second part of the uncertainty was due to residual disagreements in 
high energy data between real and simulated distributions used for calibration (correlated 
between different years). The largest contribution to the detector calibration uncertainty 
was the calibration of the HPC response. 

The uncertainty of the residual background level consisted of three parts: the precision 
of the event generators (correlated); the statistics of the simulated background events 
(uncorrelated); and the uncertainty of the background level normalisation procedure. 
The normalisation uncertainty was determined from the statistics of real data events used 
for the normalisation (uncorrelated) and from residual data-simulation disagreements in 
the distributions of variables sensitive to particular background types (correlated). The 
uncertainties of the background from different sources are presented in Table 13. 

Other (smaller) sources of systematic uncertainties were: statistics of simulated sig¬ 
nal events; QED uncertainties computed by changing the order of perturbation theory to 
which QED corrections are included in the KK generator; subtraction of r pairs accompa¬ 
nied by light pair production; reconstruction of estimation of the feed-through event 
fraction; beam energy spread; trigger efficiency; world averages of r-decay branching 
ratios; and r-polarisation at LEP II energies. 

An uncertainty specific to the measurements of the forward-backward asymmetries and 
differential cross-sections originated from the estimation of the charge misidentihcation 
probability. The systematic errors on the total non-radiative cross-section were assigned 
also to the differential cross-section. A correlation between the cos 9 bins was assumed 
for all sources except the statistics of simulated signal and simulated background. The 
charge misidentihcation uncertainty was assumed to be anti-correlated between the bins 
with cos 6* < 0 and cos^ > 0. 

A relatively small additional systematic uncertainty for the period of the broken TPC 
sector came from reconstruction, trigger efficiency and charge identification. It was 
taken into account in the procedure of averaging of the two periods of 2000. 

Cross-check with Z° data 

For an additional cross-check of the data quality we have used the data taken in 1997- 
2000 during the LEP runs near the peak of the Z resonance. After the usual run selection 
the total integrated luminosity was about 11 pb“^. Tan pairs were selected from this 
sample using a set of cuts very similar to the ones used for the high-energy analysis, and 
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the total cross-section and forward-backward asymmetry were calculated using the same 
technique as for the high energies. 

For the combined 1997-2000 sample the statistical uncertainties were Aa/a = ±1.2% 
for the cross-section and AApB = ±0.012 for the asymmetry estimation. Both values are 
similar to, or below, the systematic errors of the measurements at high energy. 

The measured value of the forward-backward asymmetry was in good agreement with 
the Standard Model expectation, while the cross-section was found to be 1.1 standard 
deviations below the expected value (only the statistical error was taken into account). 
The year-by-year studies also did not show any signihcant departure from the expecta¬ 
tions. These cross-checks give additional conhdence that the measurements performed at 
the high energies are valid within the quoted uncertainties. 

5.4 qq final states 

The analysis of the e’''e“ —qq( 7 ) process followed closely the procedure applied pre¬ 
viously to the data collected near 183 and 189 GeV and published in [4]. It benehted 
however from several improvements. In particular, more accurate generators were used 
for the simulation of the signal and background hnal states, and the enlarged Monte 
Carlo data sample allowed an improved understanding and control of residual differences 
between real and simulated distributions. 

These improvements were implemented in the analysis of the latest data sets (collected 
above 190 GeV) and stimulated a reanalysis of the published data collected near 183 
and 189 GeV. The magnitude of the changes as a result of these improvements did not 
justify reanalysing the data collected at 130-172 GeV [3], which have larger statistical 
uncertainties. 

Two cross-sections were computed at each collision energy: an inclusive cross-section 
corresponding to reduced centre-of-mass energies larger than 10% of the collision en¬ 
ergy, and a non-radiative cross-section corresponding to reduced centre-of-mass energies 
exceeding 85% of the collision energy. 

5.4.1 Analysis 
Run selection 

The cross-section computations are based on data samples collected during running 
periods where the subdetectors of prime importance for the event selection were close 
enough to their nominal operating conditions. A small fraction of the data collected 
was discarded for this reason, mainly due to inefficiencies of the TPG, of the forward 
chambers (FGA-FGB) or of the FEMG. They correspond to about 1.4% of the total 
integrated luminosity. 

Track selection 

The selection of qq events relied mainly on charged particles for which the tracks 
were required to originate from the vertex (i.e. their transverse and longitudinal impact 
parameters with respect to the nominal interaction point had to be below 4 cm and 10 cm 
respectively), to have a momentum greater than 400 MeV/c measured with a resolution 
better than 100%, and to have a reconstructed track length of at least 30 cm. For a 
small fraction of particle trajectories which were only reconstructed in the VD and ID, 
the length was required to exceed only 10 cm. 
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Figure 7: Charged multiplicity (left) and transverse energy (right) distributions for 
> qq events. The arrows indicate the cut values applied. The distributions combine 
data collected from 1997 to 2000. 











42 


DELPHI 


DELPHI 




60 70 80 90 100 110 120 130 140 150 160 


Vs’, GeV 


Figure 8: Reconstructed \/^/y/s (left) and (right) distributions for e —>• qq events. 
Data are from 1997-2000. 

Event selection 

The event selection was mainly oriented towards an efficient rejection of the back¬ 
grounds due to low multiplicity events and two-photon collisions. For this purpose, final 
states were accepted if they contained at least 7 selected tracks and if the total energy 
carried by the selected tracks exceeded 10% of the collision energy. In order to suppress 
the contamination by two-photon collisions further, the event total transverse energy was 
computed, based on the transverse momentum of each selected charged particle and on 
the electromagnetic showers reconstructed in the HPC and the FEMC with a shower 
energy above 500 MeV. Events were rejected if their total transverse energy was below 
20% of the collision energy. Figure 7 shows the distributions of charged multiplicity and 
transverse energy for candidate e+e“ —> qq events combined over the years 1997 to 2000. 
Finally, a large fraction of the residual Bhabha events was discarded by requiring the 
quadratic sum of the total energy (i.e. Erad = \J+ Eg) seen in each end-cap of the 
FEMC (computed with showers of at least 1 GeV) to be less than 90% of the beam energy. 
The sample of events satisfying the whole set of selection criteria described previously 
was thus mainly contaminated by four-fermion hnal states. 

The selection of non-radiative hnal states included the additional requirement that the 
reduced centre-of-mass energy was larger than 85% of the collision energy (i.e. a/s' > 
0.85 Vs). 

y/~P determination 

The determination of the reduced energy of each event (VV) was mainly based on 
hadronic jets. A hrst step consisted in reconstructing jets according to the DURHAM 
clustering algorithm [30], using the selected charged particles as well as the neutral par¬ 
ticles with energy larger than 1.5 GeV. A constrained ht was then performed, imposing 
energy and momentum conservation and assuming that an ISR photon was emitted along 
the beam lines. The free parameters of the ht were the jet directions and the ISR photon 
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^e“ ^ 

qq 






Energy point (GeV) 

183 

189 

192 

196 

200 

202 

205 

207 

Energy (GeV) 

182.7 

188.6 

191.6 

195.5 

199.5 

201.6 

204.9 

206.5 

Lumi (pb“^) 

53.5 

155.0 

25.1 

76.1 

83.0 

40.3 

81.9 

136.9 

No. Events 

5859 

15582 

2433 

7241 

7452 

3573 

6792 

10982 

Efficiency (%) 

92.3 

92.3 

92.4 

92.3 

92.5 

92.4 

92.1 

92.2 

Background (%) 

26.7 

30.3 

31.2 

32.3 

34.4 

35.2 

36.6 

37.6 


Table 18: The table shows the actual centre-of-mass energy and luminosity analysed at 
each energy point, the number of events selected in the inclusive samples and the efficien¬ 
cies and backgrounds for the non-radiative selection. 

energy (E^). The value of was then derived from the htted value of according to 
the following expression: s' = s — 2E^y/s. 

The quality of the agreement between the real and simulated distributions of y/^, 
which is essential for an accurate determination of the non-radiative event selection effi¬ 
ciency, was reasonable, as illustrated in Figure 8. The differences observed for values 
close to ^/s and to the Z mass reflect residual discrepancies between the real and sim¬ 
ulated yf^ resolutions. They have minor consequences on the analysis, which is mainly 
sensitive to differences occurring near the cut value of 0.85 ^/s. 

Estimation of selection efficiency and backgrounds 

The selection efficiency was derived at each collision energy from the simulated sample 
of signal events satisfying the criteria above. The sample was generated with the KK 
program [17], the event hadronisation being performed with the PYTHIA algorithm [19, 
20]. The selection efficiency values obtained for the inclusive sample varied from about 
84% to 78% for increasing values of the collision energy. It remained near 92% at all 
collision energies for the non-radiative sample. Values of the selection efficiency at various 
centre-of-mass energies may be found in Table 18. The fraction of events generated below 
85% of ^/s and reconstructed above the cut value represented less than 10% of the non- 
radiative sample. 

The residual background events contaminating the selected non-radiative {inclusive) 
sample originated mainly from WW pairs, which amounted to more than 90% (70%) of 
the total background. The inclusive sample contained also significant contributions from 
e’''e“ —ZZ, e^e“ —e’''e“ff and two-photon collision events. Other backgrounds, like 
those containing lepton pairs of opposite charge, amounted only to a few percent of the 
total residual background. 

In order to assess the magnitude of the residual contamination of the non-radiative 
sample by WW events, which accounted for up to 1/3 of the sample size, a dedicated study 
of the accuracy of the WW generator program (WPHACT) was performed. Variables 
reflecting the event shapes were multiplied to create a multi-variable selection parameter 
(see Figure 9), which was applied to the data in order to reject a substantial fraction of the 
background. The variables included: the energies of the most and least energetic jets; the 
minimal interjet angle; the broadening of the narrowest jet; the value of the clusterization 
parameter of the LUCLUS [20] algorithm at which the event changed from 3 to 4 jets. 
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Figure 9: Distribution of the multi-variable selection parameter used to tag WW events 
in the non-radiative sample. The light area represents the signal simulation, the darker 
area stands for the WW —^ qqqq background, and the darkest area displays the remaining 
backgrounds. The distributions combine data collected from 1997 to 2000. 

About 60-70% of the residual WW events were discarded in this way, while the signal 
selection efficiency dropped by a few percent. The cross-sections computed before and 
after applying these additional rejection criteria differed by only small amounts (0.2-0.3 
pb), fully compatible with the statistical uncertainty affecting the comparisons. 


Forward region simulation 

The simulated rate of tracks belonging to multi-hadronic final states and reconstructed 
in the forward regions of the detector (i.e. with a polar angle between 9° and 30° or 
between 150° and 171°) underestimated the observed rate by up to 10%. This deficit was 
regarded as a consequence of the approximations made in the description of the forward 
material of the detector at the simulation level. 

In order to account for the deficit, an artificial inefficiency was applied to the real tracks 
reconstructed in the forward regions. It consisted in discarding randomly some of these 
tracks, with weights extracted from the ratio between the real and simulated track polar 
angle distributions in four different momentum bins. These corrections were extracted 
from the calibration data collected each year near the Z resonance peak and applied to the 
high-energy data of the same year. The corresponding change of the inclusive and non- 
radiative cross-sections was —1.2% on average. The associated systematic uncertainty 
was estimated as one third of the change (i.e. 0.4%). 
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An alternative treatment was used as a cross-check. It consisted in assigning weights 
to the simulated forward tracks in order to reproduce the excess of the real ones. This ap¬ 
proach resulted in values of the inclusive cross-section well consistent with those obtained 
with the previous method, which was however preferred for the cross-section determina¬ 
tion as it was signihcantly simpler to implement in the reconstruction algorithm. 

Impact of the TPC instability in the year 2000 

The failure of one of the TPC sectors in the year 2000 translated into a small loss in 
the reconstruction efficiency for tracks traversing the sector affected, which was estimated 
to occur in about 40% of multihadronic events. According to the event simulation, 
the corresponding drop of the event selection efficiency for the inclusive sample was 
moderate (—0.6%). Several cross-checks were performed in order to test how accurately 
the simulation reproduced the modihed track and jet reconstructions. While the rate of 
real tracks was reproduced quite well, residual differences were found between real and 
simulated jet parameters, which were accounted for by a specific systematic uncertainty. 
The latter was estimated as ±0.1% for the inclusive sample selection efficiency and ±0.5% 
for the non-radiative sample selection efficiency. 

Subtraction of interference between ISR and FSR 

Computations performed with the ZFITTER package showed that the inclusion of the 
ISR-FSR interference in the computation of the non-radiative cross-section diminishes 
its magnitude by almost 0.5%. 

Since the KK generator used for the simulation of multi-hadronic hnal states did not 
account for interference between initial and hnal state photon radiation, and since the 
expected change in the cross-section is not uniform, peaking at | cos6*| = 1, where the 
efficiency is smallest, the quark polar angle distribution was reweighted by the relative 
difference in the differential cross-sections computed from ZFITTER, with and without 
ISR-FSR interference. The change in the selection efficiency was computed, and found 
to be negligible. Thus the measured non-radiative cross-sections were corrected by the 
full 0.5% correction. 


5.4.2 Results 


Determination of the cross-sections 


The e+e —> qq cross-sections were derived from the expression: 


_ {Nsel 

(^qq — 


Nbq) (1 - /) 

e C 


(14) 


where Ngei stands for the number of events selected, Wg is the number of background 
events expected, e is the selection efficiency and C stands for the integrated luminosity. 
The term (1 — /) applies only to the non-radiative cross-section, the parameter / ex¬ 
pressing the feed-up by events produced at a reduced centre-of-mass energy below 85% 
of 

The cross-section values found at each energy point are presented in Table 19 and 
shown in Figure 10, together with the Standard Model predictions computed with the 
ZFITTER program. The non-radiative cross-sections at energies from 130 GeV to 
172 GeV supersede those of [3] because of the correction for ISR-FSR interference men¬ 
tioned in Section 5.4.1. 
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e^e ^ qq 


Vs'/y/s>Q.lQ 

W/yTs > 0.85 


W/y/s > 0.10 

y/V/y/s > 0.85 

(GeV) 

CJqq (pb) 

CJqq (pb) 

(GeV) 

CJqq (pb) 

CJqq (pb) 

130 

328.4 ±11.3 ±3.7 

82.4 ±5.2 ±2.6 

192 

92.9 ±2.4 ±1.0 

22.1 ± 1.1 ± 0.3 

(328.1) 

(82.5) 

(93.4) 

(21.2) 

136 

259.6 ±10.0 ±3.1 

65.3 ±4.7 ±2.1 

196 

91.1 ± 1.4 ±0.9 

21.2 ±0.6 ±0.3 

(272.0) 

(66.4) 

(88.7) 

(20.1) 

161 

158.3 ±4.4 ±2.0 

41.0 ± 2.1 ± 1.3 

200 

85.2 ±1.3 ±0.9 

19.5 ±0.6 ±0.3 

(151.4) 

(35.1) 

(84.2) 

(19.0) 

172 

125.5 ±4.2 ± 1.9 

30.4 ± 1.9 ± 1.0 

202 

84.2 ±1.9 ±0.9 

18.9 ±0.8 ±0.3 

(125.1) 

(28.7) 

(82.0) 

(18.5) 

183 

107.6 ± 1.7 ± 1.0 

25.5 ±0.8 ±0.3 

205 

77.8 ±1.3 ±0.9 

17.7 ±0.6 ±0.3 

(106.0) 

(24.2) 

(78.9) 

(17.8) 

189 

96.9 ± 1.0 ±0.9 

22.6 ±0.5 ±0.3 

207 

74.7 ±1.0 ±0.8 

17.0 ±0.4 ±0.3 

(97.3) 

(22.1) 

(77.3) 

(17.4) 


Table 19: Measured values of the inelusive and non-radiative eross-seetions for the proeess 
qq at eollision energies ranging from 130 to 207 GeV. The uneertainties include 
the statistical and all systematic contributions. Values in parentheses are the Standard 
Model predictions computed with the ZFITTER program, and are estimated to have a 
precision o/±0.26%. 
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Systematic uncertainty 

Most of the systematic uncertainties on the selection efficiency resulted from small 
residual discrepancies between real and simulated distributions. Because of these differ¬ 
ences, the selection cuts were expected to have slightly different effects on the simulated 
distributions than on the real ones. The associated uncertainties were estimated from 
the observed changes of the cross-section when varying the cut values. They added up to 
a total contribution of 1.0%. 

The systematic uncertainty related to the ISR modelling was estimated by comparing 
the selection efficiencies extracted from two different simulation samples, produced either 
with the KK or with the PYTHIA generators. The observed difference in the selection 
efficiency was very small, i.e. 0.1%, and was converted into a systematic nncertainty of 
identical size. 

The uncertainty related to the accnracy of the fragmentation modelling in the Monte 
Carlo was derived from the comparison between a sample where the hadronisation was 
performed with PYTHIA and a sample where it was done with ARIADNE (version 
4.08) [31] , with parameters tnned to DELPHI data [21]. The difference between the two 
valnes of the selection efficiency (0.1-0.25%) was taken as the corresponding systematic 
uncertainty. 

A 0.7% systematic uncertainty was assigned to the non-radiative cross-section as a 
conseqnence of the cnt on Rs value follows from the resolution attributed to the 
reconstructed value of 

The systematic uncertainty on the residual backgrounds accounts for the theoretical 
precision of the generators used and for the accuracy on the selection efficiency associ¬ 
ated to each hnal state. The uncertainty affecting the non-radiative cross-section was 
completely dominated by the precision on the WW backgronnd estimation, while the 
inclusive cross-section was also affected by the nncertainties related to the snbtraction 
of residnal contaminations dne to e+e“ —ZZ, e’''e“ —> e’''e“ff and two-photon collision 
hnal states. 

The quadratic combination of the contribntions to the systematic nncertainties de¬ 
scribed above translate into a total nncertainty on the selection efficiency amonnting to 
1 .1% for the inclusive cross-section and 1.3% for the non-radiative cross-sections. As for 
the residnal backgronnds, a 0.5 pb (0.2 pb) total nncertainty was assigned to the inclusive 
{non-radiative) cross-section. The correlation between systematic uncertainties assigned 
to cross-sections measured at different energy points was estimated as 80% when coming 
from the event selection, and as 100% when originating from the residual background 
subtraction. 

The breakdown of the systematic uncertainties on the e’''e“ —qq cross-sections is 
provided in Table 20. 

5.5 Summary of results 

The measurements of cross-sections and forward-backward asymmetries for non- 
radiative samples of events are compared to theoretical predictions in Figures 10 and 11. 
In both cases the measnrements are compared directly to the predictions and the ratios 
of the measurements to the prediction are shown, 

^ “ n ’ 

^SM 
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e^e - 

qq 


Source 

Aaja 

(non-rad.) 

Aa/a 

(Indus.) 

Forward corrections 

40 

40 

Selection cuts 

100 

100 

ISR modelling 

10 

10 

Fragmentation 

25 

10 

\/p cut 

70 

- 

Backgrounds 

120 

70 

Luminosity 

55 

55 

TPC sector instability 

50 

10 

Total uncertainty 

186 

140 

Correlated 

169 

124 

Uncorrelated 

79 

65 


Table 20: Breakdown of the systematie uneertainties on the qq eross-seetions for 

data taken at ^/s ~ 207 GeV. All numbers in units of 10“"^. The total uneertainties do 
not inelude the uneertainty due to the TPC instability, which applies only to a part of the 
data eolleeted in 2000. The correlated error component includes errors correlated between 
energies and channels and those correlated with other LEP experiments. 
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Figure 10: Comparisons of the measurements of cross-sections to predictions of the Stan¬ 
dard Model from BHWIDE and ZFITTER for non-radiative samples of , 

r+r“ and qq final states. The measurements are compared to the predictions shown as 
curves (left) and as values of 71 the ratio of the measurement to the predictions (right). 
Eor clarity the results for pCpi~ and final states are slightly displaced from the ac¬ 
tual centre-of-mass energy of the data. The errors include the statistical and experimental 
systematic errors. 
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Figure 11: Comparisons of the measurements of forward-backward asymmetries to pre¬ 
dictions of the Standard Model from BHWIDE and ZFITTER for non-radiative samples 
ofe~^e~, and final states. The measurements are compared to the predictions 

shown as curves (left) and as values of IZ the ratio of the measurement to the predic¬ 
tions (right). Eor clarity the results for pCpi~ and final states are slightly displaced 
from the actual centre-of-mass energy of the data. The errors include the statistical and 
experimental systematic errors. 





51 


Meas. 

Final 

States 

(77) ± (stat) ± (syst) ± (theory) 

a 

e’^e” 

1.0006 ± 0.0086 ± 0.0077 ± 0.0200 


p^p~ 

0.9961 ± 0.0244 ± 0.0062 ± 0.0040 


r+r“ 

0.9852 ± 0.0341 ± 0.0203 ± 0.0040 


qq 

1.0256 ± 0.0103 ± 0.0130 ± 0.0026 


1+1- 

0.9930 ± 0.0200 ± 0.0074 ± 0.0040 

Afb 

e+e” 

0.9896 ± 0.0061 ± 0.0043 ± 0.0244 


/i+p- 

1.0120 ± 0.0335 ± 0.0010 ± 0.0068 



1.1121 ± 0.0405 ± 0.0156 ± 0.0068 


1+1- 

1.0494 ± 0.0259 ± 0.0059 ± 0.0068 


Table 21: Average values of TZ for eross-sections and forward-baekward asymmetries for 
, r+r“ and qq final states. The errors on the averages, (71), are the statistieal, 
experimental systematie and theoretieal uneertainties. 

where O is the measured observable and Osm is the prediction for the observable. A 
useful guide to the compatibility of the different sets of measurements with the theoretical 
predictions is to average the ratios of the measured results to the predictions over all 
centre-of-mass energies, {TV). Such averages have been calculated for the measurements 
made on non-radiative samples of events: 

• cross-sections for e+e“, joVjjT, qq hnal states and 1^1“, an average of the 

combination of the fi~and hnal states; and 

• forward-backward asymmetries for e^e“, pApT, and I’^l” hnal states; and 

• bin-by-bin averages of the diherential cross-sections for p~ and hnal states. 

The averages were made using the BLUE technique [29] so that correlations between 
systematic errors in the diherent averages were taken into account. The combination 
method makes it easy to identify the statistical and systematic uncertainties on {TV), 
for each set of measurements. For the cross-section and forward-backward asymmetry 
measurements the averages are made over all centre-of-mass energies from 130 to 207 GeV. 
Data taken at lower energies, from 130 to 172 GeV, are taken from [3], except for the qq 
hnal state where the updated values of Table 19 are used. The experimental systematic 
errors for the lower energy data, with the exception of luminosity errors, are taken to be 
correlated amongst themselves for each measurement, and uncorrelated with the errors 
at higher energies. Theoretical uncertainties on the cross-sections and asymmetries are 
taken from [32]. For the diherential cross-section the averages are made over energies 
from 183-207 GeV, and for the p'^p~ and hnal states the statistical error is taken 
from the expected error. The relative theoretical uncertainties on the diherential cross- 
sections are taken to be the same as the relative uncertainties on the cross-sections as 
given in [32], and assumed to be fully correlated between bins. 

Results for the averages of (77) for cross-sections and forward-backward asymmetries 
are given in Table 21. They indicate satisfactory agreement between the data and the 
predictions, with the largest deviation from the expectation being approximately 2.6 stan¬ 
dard deviations. In most cases the uncertainties on the averages are dominated by statisti¬ 
cal errors, except for the cross-section and forward-backward asymmetry for e’''e“ —> e’''e“ 
where the theoretical errors dominate, and for the cross-section for e+e“ —> qq where the 
experimental systematic error dominates. 
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e+( 

a 

^ e+< 





cos 9 

(77)± (stat)±(syst)±(theory) 

[-0.72,-0.54] 

1.106 

± 

0.091 

± 

0.049 

± 

0.020 

[-0.54,-0.36] 

1.077 

± 

0.074 

± 

0.014 

± 

0.020 

[-0.36,-0.18] 

0.951 

± 

0.060 

± 

0.008 

± 

0.020 

[-0.18, 0.00] 

1.113 

± 

0.049 

± 

0.026 

± 

0.020 

[ 0.00, 0.09] 

1.005 

± 

0.047 

± 

0.027 

± 

0.020 

[ 0.09, 0.18] 

0.999 

± 

0.055 

± 

0.007 

± 

0.020 

[ 0.18, 0.27] 

0.958 

± 

0.046 

± 

0.007 

± 

0.020 

[ 0.27, 0.36] 

0.943 

± 

0.038 

± 

0.007 

± 

0.020 

[ 0.36, 0.45] 

1.030 

± 

0.033 

± 

0.008 

± 

0.020 

[ 0.45, 0.54] 

0.974 

± 

0.026 

± 

0.008 

± 

0.020 

[ 0.54, 0.63] 

1.007 

± 

0.021 

± 

0.008 

± 

0.020 

[ 0.63, 0.72] 

0.999 

± 

0.016 

± 

0.008 

± 

0.020 





—/i+/i 




cos 6* 

(77)± (stat)±(syst)±(theory) 

[-0.97,-0.80] 

1.166 

± 

0.150 

± 

0.009 

± 

0.004 

[-0.80,-0.60] 

0.964 

± 

0.138 

± 

0.008 

± 

0.004 

[-0.60,-0.40] 

1.092 

± 

0.127 

± 

0.008 

± 

0.004 

[-0.40,-0.20] 

0.761 

± 

0.112 

± 

0.007 

± 

0.004 

[-0.20, 0.00] 

1.031 

± 

0.098 

± 

0.007 

± 

0.004 

[ 0.00, 0.20] 

1.004 

± 

0.085 

± 

0.007 

± 

0.004 

[ 0.20, 0.40] 

1.067 

± 

0.072 

± 

0.007 

± 

0.004 

[ 0.40, 0.60] 

1.083 

± 

0.063 

± 

0.006 

± 

0.004 

[ 0.60, 0.80] 

0.933 

± 

0.056 

± 

0.006 

± 

0.004 

[ 0.80, 0.97] 

1.003 

± 

0.054 

± 

0.006 

± 

0.004 





—> r+r 




cos 6* 

(77)± (stat)±(syst)±(theory) 

[-0.96,-0.80] 

0.397 

± 

0.291 

± 

0.025 

± 

0.004 

[-0.80,-0.60] 

0.757 

± 

0.221 

± 

0.031 

± 

0.004 

[-0.60,-0.40] 

0.970 

± 

0.166 

± 

0.023 

± 

0.004 

[-0.40,-0.20] 

0.898 

± 

0.144 

± 

0.022 

± 

0.004 

[-0.20, 0.00] 

0.972 

± 

0.133 

± 

0.023 

± 

0.004 

[ 0.00, 0.20] 

1.122 

± 

0.111 

± 

0.021 

± 

0.004 

[ 0.20, 0.40] 

1.106 

± 

0.089 

± 

0.020 

± 

0.004 

[ 0.40, 0.60] 

0.967 

± 

0.078 

± 

0.021 

± 

0.004 

[ 0.60, 0.80] 

0.875 

± 

0.086 

± 

0.022 

± 

0.004 

[ 0.80, 0.96] 

1.027 

± 

0.098 

± 

0.021 

± 

0.004 


Table 22: Bin-by-bin average values of 71 for differential cross-sections for 

and T^T~ final states. The errors on the averages, (71), are the statistical, experimental 

systematic and theoretical uncertainties. 
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DELPHI 




-1 -0.5 0 0.5 1 


-1 -0.5 0 0.5 1 


COS0 


COS0 


Figure 12: Bin-by-bin averages, (71), for jj,'^iJi~ and final states. The lower plots 
show the values of {71) while the upper plots show the SM da/dcos6 multiplied by {7Z) at 
the luminosity weighted centre-of-mass energy for each measurement channel. 
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Results for the bin-by-bin averages of ratios of the measured differential cross-sections 
to the predictions are given in Table 22. These are shown in Figure 12, in which the aver¬ 
aged ratios have been expressed as a differential cross-section at the luminosity weighted 
centre-of-mass energy for the and hnal states by multiplying the averaged 

ratios, (TZ), by the predicted differential cross-sections in each bin at the appropriate 
centre-of-mass energy. The results show good agreement with the expectation. 

Overall the data show good agreement between the measurements and the predic¬ 
tions of the Standard Model from ZFITTER and BHWIDE. In the following section the 
measurements are interpreted in a variety of models which allow for physics beyond the 
Standard Model. Since each of these models predicts specihc behaviours for either the 
differential cross-section, or alternatively cross-sections and forward-backward asymme¬ 
tries, as a function of centre-of-mass energy, the predictions of the models are compared 
directly to the individual measurements rather than to the averaged (TZ) values given in 
this section. 


6 Interpretation 

The results of the measurements of cross-sections, forward-backward asymmetries and 
angular distributions presented in Section 5 have been interpreted in a variety of scenarios. 
The data were used to determine the parameters of the S-matrix formalism for e’''e“ —> ff, 
as well as to investigate a variety of models which include explicit forms of physics beyond 
the Standard Model^: models with Z bosons, contact interactions, models which include 
the exchange of gravitons in large extra dimensions and models which consider possible s 
or t channel sneutrino exchange in R-parity violating supersymmetry. For the S-matrix 
ansatz and the search for Z bosons the LEP II data reported here were used together 
with measurements from LEP I. For other studies, the LEP II data alone were used. 

In the following subsections, each of these interpretations is discussed. In all cases the 
theoretical basis of the model is summarised and the relevant parameters are identihed. 
The methods of fitting parameters of the models to the data are discussed. The theoretical 
uncertainties on the Standard Model predictions for cross-sections and asymmetries, as 
given in Section 5, are explicitly included in the fits. In most cases, the relevant papers 
discuss the predictions of the models only at Born level, however, to confront the models 
with data it is necessary to take into account QED radiative corrections. The methods 
used to implement radiative corrections are described for each model. In the absence 
of evidence for physics beyond the Standard Model, limits on the physical parameters 
of the models can be derived. There is no universally agreed way to compute limits, 
particularly when measurements he close to a physical boundary. Where possible results 
of the interpretations are first quoted as measurements, and the methods used to derive 
limits from these results are discussed. 

In principle the introduction of processes beyond the Standard Model will affect the 
cross-section for the channel e^e“ —e^e“ at small angles, on which the measurement of 
the luminosity is based. However for values of the parameters of the models discussed in 
this section, the changes to the cross-section in this region are negligible compared with 
the experimental and theoretical uncertainties discussed in Section 4. 

"^In the comparison of the data with the models, it has been assumed that the one-loop electroweak corrections to the 
observables are those given by the Standard Model. It has been pointed out [33] that this may not always be the case. 
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6.1 S-matrix ansatz 


6.1.1 Theory 

The S-matrix formalism [34] is a rigorous semi model-independent approach to describe 
the cross-sections and the forward-backward asymmetries in e’''e“ annihilations. In this 
model, the cross-sections can be parametrised as follows: 




g|; j^(g - 

/ (s-MzV)2 + MzVrz'_ 


with 


a = tot, fb 
f = had, e, , 


(15) 


while the forward-backward asymmetries are given by: 



(16) 


The mass Mz and width Tz of the Z in the S-matrix formalism are related to the values 
obtained from Standard Model hts, Mz and Tz, in which an s dependent width term is 


included by 

Mz = Mz-\/l tTzVMz^ ~ Mz+34.20 MeV/c^ 
Tz = Fz^l + rz'/Mz' ^ rz+0.94MeV. 


(17) 


The parameters if and jf respectively scale the Z exchange and the yZ interference 
contributions to the total cross-section and forward-backward asymmetries; they are func¬ 
tions of the effective vector and axial-vector couplings of the fermions. The contribution 
gf of the pure 7 exchange was hxed to the value predicted by QED in all hts. The 
photonic virtual and bremsstrahlung corrections are included through the convolution of 
eqn. 15 with the photonic hux function. 


6.1.2 Results 

Published measurements from LEP I [1,2] and the runs above the Z in 1995-97 [3] 
and the results described in Section 5 were analysed in the framework of the S-matrix 
approach, achieving a substantial improvement in the precision of the 7 Z interference 
compared to the accuracy obtained from the Z data alone [1,2], and updating the results 
presented in reference [3]. 

Fits to the hadronic and leptonic cross-sections and leptonic forward-backward asym¬ 
metries were carried out in this framework using the corresponding branch of the 
ZFITTER/SMATASY 6.36 [26,34,35] program. Data on e’''e“ —e’'“e“ at LEP II are 
not used in the ht due to the large f-channel contribution to the measurements, which is 
not described by the S-matrix formalism. In the fits, LEP I and LEP II measurements 
were assumed to be uncorrelated. Results for the mass and the width are quoted in terms 
of Mz and Pz. 

The results of the hts are presented in Table 23. Using the LEP I data only, the 
amounted to 162.0 {ndof = 161) for the 16-parameter ht (i.e. without assuming lepton 
universality), and to 176.1 {ndof = 169) for the 8 -parameter ht (where lepton universality 
was assumed), the number of htted points being 177. Using the combined LEP I and 
LEP II data, the amounted to 245.1 {ndof = 221) and 256.1 {ndof = 229) for the 
16-parameter ht and the 8 -parameter ht, respectively, the number of htted points being 
237. The correlation coefficients between the free parameters of the 16- and 8 -parameter 
hts for the LEP I and LEP I -|- LEP II data are shown in Tables 24, 25, 26 and 27. 
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The data support the hypothesis of lepton universality. Overall, the measurements are 
in good agreement with the Standard Model predictions. 

The correlations between the parameters Mz and j};° ^ for 8-parameter hts to LEP I 
and LEP I -|- LEP II data are shown in Figure 13. It can be seen that a signihcant 
improvement on the precision on the hadronic interference parameter, is obtained 
when the high energy data are included in the £t. The fitted value of Mz is consistent 
with the value obtained from a 5 parameter. Standard Model, fit to LEP I data [2], 
Mz = 91.1863±0.0028 GeV/c^, and the value of is consistent with, though somewhat 
above, the Standard Model expectation = 0.2201]';o;oo3g, where the uncertainties 
correspond to uncertainties in the input parameters. 

6.2 Z' Bosons 

6.2.1 Theory 

Many theories which are more general than the Standard Model predict the existence 
of additional heavy gauge bosons. The consequences of several of these models were 
investigated, complemented by a model-independent fit to the leptonic data. 

The existence of an additional heavy gauge boson Z can be parametrised by the 
mass of the boson M^' and by its couplings to fermions. In addition, a possible mixing 
between Z' and the standard Z, represented by a mixing angle ©zz'’ taken into 

account [36,37]. In order to deal with a restricted number of free parameters, it is useful 
to consider specific Z models with well-defined couplings. Popular models are: 

• The Eq model [38,39] is based on a symmetry breaking of the Eq GUT where two 
gauge groups U(l)^ and U{1)^ are introduced, 

^ 2 ' = cos ©6sin ©6 . (18) 

The free parameter, ©g, of this model is the mixing of the y and ip fields to form 
the Z'. Usual choices of ©e are ©e = 0, n/2, — arctan ^J5/3 ( y, ip and rj model); 

• The L-R model [40,41] includes a right-handed S'U(2)jj extension to the Standard 
Model gauge group SU{2)l ® U{1). The free parameter aiR describes the coupling 
of the heavy bosons to fermions. aiR can take values between ^2/3 < aiR < 

\J cot^ 9w — 1 (~ 1.53), where 9w is the weak mixing angle. 

In a more general approach, the Z! boson is directly described in terms of its couplings 
a'j and v'j [42]. Off the Z resonance, pair production is only sensitive to the normalised 
couplings and . As a consequence, the couplings and the mass of the Z! boson 
cannot be measured independently. The normalised couplings are: 



6.2.2 Results 

The non-radiative hadronic total cross-sections and and r+r“ total cross-sections 
and asymmetries presented here together with existing data from LEP I [1,2] and LEP 
II at y/s ~ 130-172 GeV [3] were used to fit the data to models including additional Z 
bosons. 
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LEP I 

LEP I ± LEP II 

SM 

Mz 

91.1936 ±0.0112 

91.1808 ± 0.0094 

91.1844 ±0.0036 

91.1831 ± 0.0034 

- 

Tz 

2.4861 ± 0.0048 

2.4886 ± 0.0046 

2.4894 ± 0.0041 

2.4893 ± 0.0041 

2.497 

tot 

thad 

2.9490 ±0.0110 

2.9543 ±0.0107 

2.9567 ± 0.0096 

2.9564 ± 0.0095 

2.966 

tot 

0.14092 ±0.00095 


0.14129 ±0.00091 



tot 

0.14274 ±0.00072 


0.14301 ±0.00067 



tot 

t T 

0.14161 ±0.00100 


0.14204 ±0.00096 



T-tot 

Tl 


0.14230 ±0.00062 


0.14240 ±0.00058 

0.1427 

•tot 

Jhad 

-0.21 ±0.64 

0.51 ±0.55 

0.44 ±0.12 

0.47 ±0.12 

0.22 

•tot 

Je 

-0.094 ± 0.075 


-0.042 ± 0.048 



•tot 

Jm 

0.056 ± 0.042 


0.029 ±0.019 



•tot 

Jr 

0.004 ± 0.046 


-0.013 ±0.026 



•tot 

Jl 


0.047 ±0.037 


0.010 ±0.015 

0.004 


0.00306 ± 0.00091 


0.00298 ± 0.00090 




0.00275 ± 0.00051 


0.00286 ± 0.00049 



r? 

0.00406 ± 0.00072 


0.00428 ± 0.00070 





0.00304 ± 0.00038 


0.00327 ± 0.00037 

0.00273 

j? 

0.802 ± 0.075 


0.805 ± 0.075 



Jm 

0.711 ±0.037 


0.802 ± 0.024 



j? 

0.707 ± 0.047 


0.832 ± 0.031 





0.725 ± 0.027 


0.811 ±0.018 

0.799 


Table 23: Results of the 16- and 8-parameter fits to the LEP I only and combined LEP 
I -h LEP II data. Also shown are the Standard Model predictions for the fit parameters. 
Mz ( Tzj are measured in units of GeV/c^ (GeV), all other quantities are dimensionless. 
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Table 26: Correlation matrix for the 8-parameter fit to LEP I data. 



Table 27: Correlation matrix for the 8-parameter fit to LEP I + LEP II data. 
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Figure 13: Contour plot in the plane. The dashed curve shows the region 

accepted at the 68% confidence level from an 8-parameter fit to data taken at the energies 
around the Z; the solid curve shows the region accepted at the same confidence level when 
the high energy data are also included in the fit. The band labelled 5-par fit is the fit to 
Mz obtained from a 5-parameter, Standard Model, fit to LEP I data. The narrow band 
labelled SM is the Standard Model expectation for as calculated using ZFITTER. 
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Fits were made to the mass of Z\ M^.', the mass of the Z, Mz, and to the mixing 
angle between the two bosonic helds, ©zz'’ ^ different models referred to as X) t/’) V 

and L-R as described in 6.2.1. The theoretical predictions came from the ZEFIT package 
(version 6.10) [36] together with the ZFITTER program (version 6.10). The program 
ZEFIT provides predictions for the cross-sections and forward-backward asymmetries for 
each model as a function of Mz, M^, ©zzS masses of the Higgs boson, Mh, and the 

top quark, mt, the strong coupling constant, and the Z'-model parameters ©e or aiR. 

For the L-R model aiji was set to 1.1. In order to reduce the number of free parameters, 
the following input parameters were used: mt = 175 GeV/c^, Mh = 150 GeV/c^, = 
0.118. 

The correlations between the experimental errors were taken into account when the 

was calculated between the predictions and the measurements. The most important 
correlation is between the errors on the luminosity for the cross-section measurements. 
Gorrelations between LEP I and LEP II measurements are very small and therefore 
neglected. 

No evidence was found for the existence of a Z -boson in any of the models. The 
htted value of Mz was found to be in agreement with the Standard Model value. Two- 
dimensional exclusion contours at 95% conhdence level, obtained with > Xmin + 
5.99 [43] in the M 2 .'-© 2 z' plane, were made. The allowed regions for M^' and 
shown in Figure 14. The one-dimensional limits at 95% conhdence level for both M^' 
and © 22 ', obtained with > Xmin + 3.84, are shown in Table 28. We repeat that these 
limits are based on the assumption that the electroweak one-loop corrections are those 
of the SM. In particular for the L-R symmetric model, the authors of reference [33] have 
shown that inclusion of the extra gauge bosons in the corrections can affect the limits on 
the mixing angle. The limits for the Z mass range from 360 to 545 GeV/c^, an increase 
of between 50 and 125 GeV/c^ on the limits presented in [4], depending on the model. 

In addition the Sequential Standard Model [44] has been considered. This model 
proposes the existence of a Z with exactly the same coupling to fermions as the Standard 
Model Z. A limit of M^' > 1305 GeV/c^ is found at 95% conhdence level, an increase of 
595 GeV/on the limit presented in [4]. 

Mo del-independent hts were performed to the leptonic cross-sections and forward- 
backward asymmetries, for the leptonic couplings of a Z , ajf and , normalised for the 
overall coupling scale and the mass of the Z [42]. 

Several values of the mass of the Z were considered (i.e. 300, 500 and 1000 GeV/c^), 
and the ZZ -mixing was neglected. The 95% conhdence level limits on the normalised 
couplings are \aff\ < 0.19 and \vl/\ < 0.19, an increase of 0.04 and a decrease of 0.03, 
respectively, on limits given in [4]. 

6.3 Contact interactions 

6.3.1 Theory 

Gontact interactions between initial and hnal state fermionic currents provide a rather 
general description of the low energy behaviour of new physics with a characteristic high 
energy scale. Following reference [45], these interactions are parameterised by an ehective 
Lagrangian, Cgff, of the form: 

\ ’ i,j=L,R 


( 20 ) 
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Figure 14: The allowed domains in the plane for the x, 'f’, V ^-R models. 

Values outside the contours are excluded at 95% confidence level. 


Model 

X 


7] 

L-R 

(GeV/c2) 

545 

475 

360 

455 

Qhmit (mrad) 

3.1 

2.7 

9.2 

2.8 


Table 28: 95% confidence level lower limits on the T! mass and upper limits on the ZZ' 
mixing angle within the x, fi’, 9 o,nd L-R models. 
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Model 


Vrr 

Vlr 

hRL 

LL± 

±1 

0 

0 

0 

RR± 

0 

±1 

0 

0 

VV± 

±1 

±1 

±1 

±1 

AA± 

±1 

±1 

Tl 

Tl 

LR± 

0 

0 

±1 

0 

RL± 

0 

0 

0 

±1 

V0± 

±1 

±1 

0 

0 

A0=^ 

0 

0 

±1 

±1 


Table 29: Choices of rjij for different contact interaction models, LL, RR etc. 
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where j^Ti is taken to be 1 by convention, <5 = 1(0) for / = e(/ ^ e), 77 *^ = ±1 or 0, A 
is the scale of the contact interactions^, e* and fj are left or right-handed spinors. This 
effective Lagrangian is added to the Standard Model Lagrangian. By assuming different 
helicity couplings between the initial-state and hnal-state currents and either constructive 
or destructive interference with the Standard Model (according to the choice of each rjij) 
a set of different models can be dehned from this Lagrangian [46]. The values of r]ij for 
the models investigated in this paper are given in Table 29. The differential cross-section 
for scattering the outgoing fermion at an angle 6 with respect to the incident e~ direction 
is given by [47,48] 


da/d cos 6^ 


ttq;^ 


Nil 



( 21 ) 


where s, t and u are the Mandelstam variables and A"/ is the number of colours for 
fermion /. The Aij and Aij are helicity amplitudes for the scattering process [45]. When 
the helicity amplitudes are squared, 3 sets of terms arise: the first set contains purely 
Standard Model terms; the second set of terms derive from the interference between 
contact interactions and the Standard Model, these terms are proportional to 1/A^; the 
hnal set of terms are due to contact interactions alone and are proportional to 1/A^. For 
the purpose of htting contact interaction models to the data, a new parameter e = 1/A^ is 
dehned, with e = 0 in the limit that there are no contact interactions. This parameter is 
allowed to take both positive and negative values in the hts. It is worth noting that there 
is a symmetry between models with rjij = -|-1 and those with rjij = — 1. The predicted 
differential cross-section in the constructive (-I-) models is the same as the destructive (—) 
models for e~ = Therefore, starting from a model with constructive interference 

with the Standard Model, the region e > 0 represents physical values of 1/A^, while 
the region e < 0 represents physical values for the equivalent model with destructive 
interference. 


6.3.2 Results 

Cross-section and forward-backward asymmetry measurements from all LEP II centre- 
of-mass energies were compared to each of the contact interaction models mentioned 
above® using hts, considering separately the e+e“, and hnal states, and 

a simultaneous ht to all three hnal states, assuming lepton universality in the contact 
interactions to obtain curves as a function of e. The predicted tree-level results were 
corrected for QED radiation. Correlations between the systematic uncertainties were 
taken into account in the hts. Making use of the symmetry, mentioned above, between 
models with constructive and destructive interference with the Standard Model, it is 
possible to ht pairs of models by allowing e to take both positive and negative values. 

The resulting y^ as a function of e were not always Gaussian-like parabolae. The 
values of e extracted from the points with minimum y^ from each of the hts were all 
compatible with the Standard Model expectation e = 0 TeV“^, i.e. the diherences in y^ 
between the best ht points and the Standard Model point were always less than nine, 

^The choice of is somewhat arbitrary; if the coupling constant was taken to be cts much lower limits on A would be 
obtained. 

®For leptonic final states, models with only ^lr = equivalent to models with only tirl = il- 
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except in the fit for the AO model in e^e“ —> , in which case the difference in was 

found to be 11.2. 

Errors on e were derived by Ending the points above and below the best fit points for 
which the increased by nine above the minimum value, which would correspond to 
a three sigma uncertainty for Gaussian-like parabolae. These three sigma uncertainties 
were divided by three, to give one sigma errors. 95% confidence level lower limits on A, 
by integrating under the likelihood curves, C{e), were obtained from the fifs using 

C{€) = exp-^X^(e) 

over the physically allowed region of e for each choice of model^. The fitted values of e 
with upper and lower errors and the limits on A for models with constructive. A"*", and 
destructive, A“, interference with the Standard Model are given in Table 30. 


6.4 Gravity in large extra dimensions 

6.4.1 Theory 

The large difference between the electroweak scale (Mew ~ 10^ — 10^ GeV) and the 
scale at which quantum gravitational effects become strong, the Planck scale (Mpi ~ 
10^® GeV), leads to the well known “hierarchy problem”. A solution, not relying on 
supersymmetry or technicolour, that involves an effective Planck scale, Md, of (P(TeV) is 
achieved by introducing n compactified dimensions, into which spin 2 gravitons propagate, 
in addition to the 4 dimensions of standard space-time [49]. The Planck mass seen in the 
4 uncompactified dimensions, Mpi, can be expressed in terms of Mp, the effective Planck 
scale in the n + A dimensional theory, 

Mpi^ ~ R'^MD’^+^ 


where R is the size of the extra dimensions. With Mp = 1 TeV, the case where n = 1 is 
excluded as Newtonian gravitation would be modified at solar system distances, whereas 
n = 2 corresponds to a radius for extra dimensions of 0(1 mm), which is excluded 
by recent gravitational experiments which test the inverse square law of gravitational 
attraction down to (P(lOOpm) [50]. There are also severe limits from astrophysics for a 
small number (n = 2, 3) of large extra dimensions, however, higher numbers of dimensions 
are not ruled out. In addition it is possible to construct models which evade gravity and 
astrophysics bounds with a slight modification of the extra dimension scenario [51]. 

In high energy collisions at LEP and other colliders, new channels not present in the 
Standard Model would be available in which gravitons could be produced or exchanged. 
Virtual graviton exchange would affect the differential cross-section for e’''e“ —> ff, with 
the largest contributions seen at low angles with respect to the incoming electron or 
positron. Embedding the model into a string model, and identifying the effective Planck 
scale, Md, with the string scale. Mg, the differential cross-section for e’''e“ —ff with the 
inclusion of the spin 2 graviton can be expressed as [52] 


da/d cos 6^ = da/d cos 6*-f C({s, cos 6) 


X 

Mf 


-|- C2{s,cos6) 


1 2 


Mf 


'Integrating under the likelihood curve for e is equivalent to obtaining a Bayesian limit assuming a prior uniform in e 
over the physically allowed regions and zero for the unphysical regions. 
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e+e —^ e+e 

Model 


A-(TeV) 

A+(TeV) 

LL 

0 0071 

U.UU^ i_0 0135 

6.8 

5.3 

RR 

0 01 nR+0 0157 

U.UiUO_o.oi48 

6.8 

5.2 

VV 

0.00241^°! 

13.9 

11.7 

AA 

0.0176l[i:[i^i 

4.6 

4.8 

RL 

0.0035l[J:[i}i 

6.7 

5.7 

LR 

0.0035^n}i 

6.7 

5.7 

VO 

0.0038l[j:[i°^9^ 

10.1 

8.0 

AO 

0.00581^^ 

9.8 

8.5 


g —^ -J-+/j- 

Model 

el:!(TeV-2) 

A-(TeV) 

A+(TeV) 

LL 

-o.oi94;nii 

4.6 

7.9 

RR 

-0.0213l°:°l“ 

4.4 

7.6 

VV 

-0.0127l°:“i^ 

7.1 

15.8 

AA 

0.00291^°?° 

9.4 

8.8 

RL 

-0.1974^H678 

2.1 

7.9 

LR 

-0.1974^H678 

2.1 

7.9 

VO 

0 01 Q4+0-0083 
-U.UiO^_0.0069 

6.7 

11.8 

AO 

n oooq+0.0720 
-U.ZZZ^_Q 0098 

2.1 

11.6 


e+e —^ 

Model 

etlf{TeV-^) 

A-(TeV) 

A+(TeV) 

LL 

0.0019^S 

7.6 

7.3 

RR 

o.ooiotS 

7.2 

7.0 

VV 

-0.0006f^,S 

12.9 

12.2 

AA 

0.0028tHo57 

10.9 

10.1 

RL 

n oq'7'7+0.0919 
-U.ZO/ ' -0.0139 

2.0 

6.3 

LR 

n oq'7'7+0.0919 
-U.ZOi /_0.0139 

2.0 

6.3 

VO 

-O.OOlll^gll 

11.5 

10.9 

AO 

-0.23961^^66 

2.0 

9.0 


e+e“ — s- l+l~ 

Model 

el:i(TeV-2) 

A-(TeV) 

A+(TeV) 

LL 

-0.00171^°®? 

8.2 

9.1 

RR 

-0.00151°:“^? 

7.9 

8.7 

VV 

0 001 2+C' 0023 

U.UUiZ_g QQ24 

13.7 

16.5 

AA 

0 001 1 +0-0042 
U.UUii_o.o039 

12.1 

10.6 

RL 

n 0071 +0 0090 
-U.UU/^i_0 0097 

6.5 

8.7 

LR 

0 0071 +0-0090 
-U.UU/^i_0 0097 

6.5 

8.7 

VO 

0 0007+'^+*^^® 

-U.UUU 0037 

11.5 

12.6 

AO 

0 0035+*^+^^^ 
-U.UUOO_0,0048 

9.2 

12.2 


Table 30: Fitted values of e and 95% confidenee lower limits on the scale, A, of con¬ 
tact interactions in the models discussed in the text, for , 

> T~^T~ and e’''e“ ^ 1^1”, a simultaneous fit to the above, assuming lepton univer¬ 
sality in the contact interactions. The errors on e come from the statistical and systematic 
errors on the determination of the cross-sections and asymmetries. 
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while the differential cross-section for the process e’''e —e’''e can be expressed as [53] 

|2 


d(j/dcos6' = d(T/dcos6*L^ — 


A 

M* 


+ C2{s,t) 


A 

Mf 


with 0, as usual, being the polar angle of the outgoing fermion with respect to the 
direction of the incoming electron. The functions C(, C^, Cl and Cf are known [52,53]. 
For e+e“ —ff the expansion of cos 9 extends up to the fourth power. The dimensionless 
parameter A, of 0(1), is not explicitly calculable without full knowledge of the underlying 
quantum gravitational theory. It can be either positive or negative [52,54]. For the 
purposes of the fits, two cases, A = ±1, are considered. This parametrisation has no 
explicit dependence on the number of extra dimensions, n. 
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Final State 

( TeV-^) 

A 

Ms (TeV) 
[95%C.L.] 

e'''e“ 

+0.54«;g 

-1 

+1 

0.986 

0.874 

p+p- 

+0.84^1:1 

-1 

+1 

0.739 

0.695 


-3.04111 

-1 

+1 

0.578 

0.690 

1+1- 

+0.44l°:i 

-1 

+1 

0.998 

0.898 


Table 31: 95% confidence level lower limits on Ms in models of gravity in extra dimensions 
for and r+r“ final states, and for 1 ^ 1 “, a simultaneous fit to the final states 

above. 



-0.6 -0.4 -0.2 0 0.2 0.4 0.6 


COS0 

Figure 15: The deviations of the measured values of da/d cos 9 fore'^e~ —>■ e’''e“ from the 
predictions of the Standard Model, averaged over all energies. Superimposed are the pre¬ 
dicted deviations at the luminosity weighted centre-of-mass energy for gravity in large ex¬ 
tra dimensions for the 95% confidence limits. Mg =0.986 TeV for X = —1 and Mg =0.874 
TeV for A = +1. 


6.4.2 Results 

fits to the measured differential cross-sections for e’''e“, pfipT and final states 
reported in this paper for the parameter e = A/Mf were performed. In order to fit the 
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prediction to the data, the angular distributions given in [52,53] were corrected to account 
for radiative effects, dominated by initial state radiation. The corrections take the form of 
a numerical version of the radiator function contained in ZFITTER, which are computed 
assuming only Standard Model processes. Separate radiator functions were derived for 
each angular bin in several bins of V^/The tree-level new physics cross-section in a 
given bin was subsequently convoluted with the appropriate numerical radiator function 
to produce the corrected cross-section in that bin. 

For the and hnal states the statistical errors expected from the Standard 

Model were used. The correlations between the systematic errors between bins of cos 6* 
and between channels and energies were taken into account in the £t. Errors on the 
parameter e were determined in the same way as in Section 6.3.2. 

The values of e obtained are compatible with the Standard Model, i.e. e = 0 TeV”*^. 
Table 31 shows the fitted values of e and 95% confidence level lower limits on Mg. These 
limits were obtained using a method equivalent to that used to extract the limits on the 
scale. A, of contact interactions, as described in Section 6.3.2. 

The deviations of the measured angular distributions from the Standard Model pre¬ 
dictions averaged over all energies in the channel e’''e“ —> e’''e“ are shown in Figure 15 
together with the expected deviations as a function of cos 9 for the fitted best value of e 
and for the 95% lower limits on Mg for A = ±1. 


6.5 Sneutrino exchange models 

These models consider possible s- or t-channel sneutrino z/; exchange in R-parity violat¬ 
ing supersymmetry [55], which can affect the channel e’''e“ —RR. The purely leptonic 
part of the R-parity violating superpotential has the form 


njk 


LIL^lEr 


where ijk are family indices, Ll represents a left-handed leptonic superfield doublet and 
Eji corresponds to the right-handed singlet superfield of charged leptons. The coupling 
Xijk is only non-zero for combinations involving at least two generations and for i<j. 

For the channel e’''e“ —e’''e“ there are possible contributions from the s-channel pro¬ 
duction and t-channel exchange of either (A 121 % 0) or (A 131 % 0). For the channels 
e'''e“ —and e’''e“ —there is no s-channel contribution if only one of the 
XijkS is non-zero. For e’''e“ —there are t-channel contributions from either z/g 
(A 121 % 0 ), z/^ (A 122 % 0 ) or from (A 132 or A 231 % 0 ). If both A 131 % 0 and A 232 % 0 
then the s-channel production of i^r is possible. For e^e“ —> r+r“ there are t-channel 
contributions from either z/g (A 131 % 0), z/^ (A 123 or A 231 % 0) or from ( A 133 % 0). If 
both A 121 % 0 and A 233 % 0 then the s-channel production of z/^ is possible. 

All these possibilities are considered here. For a given scenario the s- or t-channel 
sneutrino exchange amplitude contribution is added to the Standard Model contribution 
as appropriate. If there is no sneutrino exchange for a specific channel then the prediction 
for that channel is just the SM value. 

In the case of s-channel sneutrino graphs, if the sneutrino mass, m~, is equal to the 
centre-of-mass energy of the e’''e“ beams, resonant sneutrino production occurs, which 
can lead to a large change in the cross-section. A lesser change in the cross-section will 
occur for m~ < due to the process of radiative return. There is some sensitivity to 
m~ just above y/s due to the finite width of the particle. It is assumed here that the 
sneutrino width is 1 GeV. 
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6.5.1 Fits to models of sneutrino exchange 

The total cross-section and forward-backward asymmetry values for the channels 
e'''e“ —e'''e“, e+e“ — jj,'^jj,~ and e’''e“ —, at each centre-of-mass energy, were used 
in the fits. The theoretical prediction consisted of Improved Born Approximation Stan¬ 
dard Model terms, plus sneutrino exchange, plus interference terms. 

All the hts considered result in values of A which are compatible with zero; so results 
are expressed as 95% conhdence limits. The hrst fits considered are to those terms which 
modify the e^e“ —> e’''e“ channel. These involve the s- and t-channel exchange of a ( 
Ai 2 i % 0) or Vr ( Ai 3 i % 0). The resulting 95% limits on A, as a function of m~, are given 
in Figure 16a. The best limits on A are obtained for the case where m~ is close to the 
actual centre-of-mass energy of the LEP collisions, but the radiative return process gives 
some sensitivity between these points. It can be seen that A greater than approximately 
0.05 can be excluded for m~ in the LEP II range of energies at the 95% conhdence level. 

Por the case that only one A value is non-zero there are only t-channel sneutrino effects 
for e^e“ —jjT and e+e“ —r^r“. The values of A obtained for the e+e“ —jjT 
channel and for the e’''e“ —channel are all consistent with zero, so results are 
expressed as 95% conhdence limits in Table 32. 

For the hts assuming that A 131 = A 232 = A, the resulting 95% conhdence limits on A, 
as a function of m~, are given in Figure 16b. A similar exclusion pattern to that obtained 
from the e’''e“ —e’''e“ channel is obtained. Values of A greater than approximately 0.07 
can be excluded for m~ in most of the LEP II range of energies at the 95% conhdence 
level. The exclusion contour for A 121 = A 233 = A is shown in Figure 16c, from which it 
can be seen that again a similar exclusion pattern is obtained. 


7 Conclusions 

Analyses of cross-sections and forward-backward asymmetries and diherential cross- 
sections in the e’''e“ —e’''e“, e’''e“ —and e’''e“ —channels, at centre-of-mass 
energies from 183 to 207 GeV, have been presented, along with measurements of cross- 
sections for e’''e“ —qq, in Sections 5.1 to 5.4. The results of the measurements are 
presented in Tables 5, 10, 15 and 19 for cross-sections and asymmetries and in Tables 6 , 
11 and 16 for diherential cross-sections. Results of the measurements of cross-sections 
and asymmetries from previous analyses [3] at centre-of-mass energies of 130 to 172 GeV 
are included in the tables. 

To compare data to the predictions of the Standard Model, averages of the ratios 
of measurements for non-radiative samples of events to the predictions of the Standard 
Model have been made over all centre-of-mass energies. The results of these averages 
are presented in Tables 21 and 22. The precision obtained on the averaged ratios for 




m~ = 100 GeV/c^ 

m~ = 200 GeV/c^ 

coupling 


^ (95% c.l.) 

^ (95% c.l.) 

A (t-chann. Ui in e+e“ - 

) 

0.19 

0.26 

A (t-chann. z/; in e’''e“ - 


0.56 

0.63 


Table 32: Limits on the various couplings A in t-channel sneutrino exchange in 

e'''e“ —and e’''e“ —r+r“. 
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DELPHI 



100 150 200 250 


m,„eutrino (GeV/c^) 



100 150 200 250 

m, neutrino (CeV/c^) 


Figure 16 : The 95 % exclusion limits for a) Ai2i(orAi3ij, as a function of m~ obtained 
from the e+e“ —>• e’''e“ channel; b) A131 = A232 = \, as a function of m~ obtained from 
the channel; c) A121 = A233 = \, as a function of m~ obtained from the 

e’''e“ —> r+r” channel. The sneutrino width is taken to be 1 GeV. 
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the cross-sections are ±2.3% for e’''e“ —i^e’''e“, dominated by theoretical uncertainties, 
±2.5% and ±4.0% for e’''e“ —and e+e“ —r+r“, respectively, dominated by sta¬ 
tistical uncertainties and ±1.7% for e^e“ —qq, dominated by experimental systematic 
uncertainties. The results are consistent with the expectations of the Standard Model. 
Full details can be found in Section 5.5. 

The measurements for non-radiative samples have been used to £t the parameters of 
the S-matrix model of e’''e“ —ff and to search for physics beyond the Standard Model 
in a number of models. In all cases correlations between the experimental uncertainties 
have been taken into account along with theoretical uncertainties. Using LEP I data [2] 
alone, the best fit to the parameter j};°j of the S-matrix Model is 0.51 ± 0.55. Including 
the LEP II measurements the best fit is 0.47 ± 0.12, consistent with, though somewhat 
above, the expectation of the Standard Model = 0.22. A complete set of results and 
description of the analysis can be found in Section 6.1. LEP I and LEP II data have 
also been used to search for Z bosons, yielding limits on the mass of such bosons which 
range from 360 to 1305 GeV/c^, depending on the model assumed. Full sets of results 
are presented in Section 6.2. 

Using measurements for non-radiative samples of leptons from the full range of centre- 
of-mass energies, limits have been derived on the scales of contact interactions and the 
scale associated with the propagation of gravitons in models with large extra dimensions; 
these are described in Sections 6.3 and 6.4. For contact interactions the limits range 
from 6.5 to 16.5 TeV, assuming lepton universality. For models of gravity in large extra 
dimensions the limits on the scales are found to be 0.898 TeV or 0.998 TeV, respectively, 
in the case when the new physics interferes either constructively or destructively with the 
Standard Model processes. Full details can be found in Section 6.4. In supersymmetric 
theories with R-parity violation, exchange of sneutrinos can affect the total and differen¬ 
tial cross-sections for e+e“ —UU. In Section 6.5 it is shown that limits at the level of 
0.05 can be set on the R-parity violating couplings for sneutrino masses in the LEP II 
range of energies. 
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